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cs| ! Abstract 
> 

lO ' We discuss the cosmological implications and the holographic dual theory of the 7D Randall- 

i> : 

■ Sundrum (RS) gravitational set-up. Adding generic matter in the bulk on the 7D gravity side, 



we study the cosmological evolution inferred by the non vanishing value of the brane-bulk energy 

exchange parameter. This analysis is achieved in detail for specific assumptions on the internal 

space evolution, including analytical considerations and numerical results. The dual theory is then 
Oh! 

D . constructed, making use of the holographic renormalization procedure. The resulting renormalized 

6D CFT is anomalous and coupled to 6D gravity plus higher order corrections. The critical point 
^ ' analysis on the brane is performed. Finally, we sketch a comparison between the two dual descrip- 

tions. We moreover generalize the AdS/CFT dual theory to the non conformal and interacting 
case, relating the energy exchange parameter of the bulk gravity description to the new interactions 
between hidden and visible sectors. 
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I. INTRODUCTION 



Many aspects of brane models have been recently developed both from the point of view 
of Standard Model building (see for example [2, 4] and references therein) and of cosmology 
([5]-[10] and references therein). In particular, Randall-Sundrum [1] cosmology [11] is among 
the most interesting insights related to brane-world cosmology (other are mirage cosmology 
[12], brane inflation [13], brane induced gravity [14]-[18], brane/anti-brane inflation [19]-[21], 
cosmologies from higher derivative corrections [22]-[24], particular examples with varying 
speed of hght [25] and cosmological evolution induced by the roUing tachion [26] , inflation in 
flux compactiflcation scenarios [27]-[29] and recent brane-world models [30]-[36]). Branes 
can be used to give origin to four dimensional gauge theories (living on a 3-brane) and hence 
localizing matter in four dimensions embedded in the ten dimensional space-time in which 
string theory lives (or eleven dimensional for M theory). On the gravity side, since gravitons 
propagate in all dimensions admitted by string theory, the validity of Newton's law seems 
to constrain the number of non compact dimensions to be four. 

Randall and Sundrum [1] proposed an alternative way of localizing gravity in four di- 
mensions without compactifying the extra dimensions. This was achieved in RSII model 
by assuming a warped extra direction, instead of a compact one (compact extra dimensions 
have been used [37] as an attempt of giving an explanation to the hierarchy problem, to 
which RSI [1] represents an alternative way out). The set-up of RSII is five dimensional 
gravity in a bulk space-time cut by a 3-brane. There is, in addition, a Z2 reflection of the 
extra dimension tranverse to the brane, having as flxed point the location of the brane. The 
result is a bound state graviton mode localized on the brane and a tower of KK modes, 
without mass gap, that give negligible corrections to the effective 4D gravity description. 
The metric solving the equations of motion for the five dimensional action is a slice of AdS^ 
copied and reflected w.r.t. the Z2 symmetry. It can be viewed as coming from the Type 
IIB string theory background for a stack of D3-branes, in the low energy effective fleld 
theory approximation, which is indeed dual to the AdS^ x near horizon geometry for a 
3-brane supergravity solution. In RS analysis only gravity on AdS^ is considered, since the 

is factored out from the anti de Sitter space, giving KK modes. The truncation of AdS^ 
with the 3-branc cuts out its boundary. 

In a recent work [38], a different string background related to RS set-up has been con- 
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sidered. The analogous analysis to the 5D RS model has been made by Bao and Lykken 
[38] in a seven dimensional anti de Sitter background rather than in the five dimensional 
original model. The background may come from the near horizon geometry of M5-branes 
in the eleven dimensional M theory, which gives AdSj x S^. As for the five dimensional 
model, the sphere is factored out and only the physics of gravity in AdSy plus KK modes is 
considered. A further step performed in [38] is to reduce AdS^ AdS^ x E^, where is a 
two dimensional internal space (namely a two-sphere or a torus) . In [38] the RS spectrum 
of KK modes gets modified and supplemental KK and winding modes appear, because of 
the E^ compactification. In [39] the 7D supergravity orbifold compactification on 3^/7,2 
is considered in the context of anomaly cancellation on the boundary of the background, 
showing that the matter contents of the theory cannot be completely generic. 

On the cosmological side, RS models can give new descriptions of the cosmological evolu- 
tion of our universe. A realistic model should be able to explain the existence of dark energy 
and the nature of dark matter, early time infiation and eventually the exit from this phase, 
as well as the late time acceleration coming from the present observations — additional 
issues are related to the cosmological constant, temperature anisotropics, etc. (see [40]-[42] 
for recent reviews on the observed cosmology). RS cosmologies have been studied and found 
to exhibit some of these features. Generally, brane-world cosmological models should take 
account of the energy exchange between brane and bulk that naturally arises because of the 
non factorized extra direction. The implications of this energy exchange has been analyzed 
in [43]. The authors propose some scenaria describing the cosmological evolution of a uni- 
verse with two accelerating phases, as we expect from experimental data. Moerover, most of 
the fixed points were shown to be stable. Earlier attempts include [44]. Subsequent papers 
[45]-[57] have been written on the subject, also finding new solutions, some of which exact. 

The aim of this paper is to investigate the 7D RS cosmology with brane-bulk energy 
exchange and to explore the model from the holographic point of view, making an explicit 
comparison between the two descriptions. Our starting point is gravity in the 7D bulk cut 
by a five-brane and with the usual RS Z2 reflection plus a generic matter term on the 
brane. In order to study the cosmological evolution of the brane-world, the ansatz for the 
metric is time dependent. Besides, the direction tranverse to the 5-brane is the warped 
direction characterizing RS models. Unlike in [38], we have different warp factors for the 
3D extended space and for the two dimensional compact internal space. It is worth noticing 



3 



that the gravitational couphng constant of the 4D space-time is dynamically related to the 
7D Newton constant, since the compact space volume is generally time dependent. Indeed, 
the 4D energy density and pressures are also dynamical functions of the density and the 
pressures defined on the brane. We analyze the generic Friedmann and (non) conservation 
equations, also including the energy exchange terms, in order to get the expressions for 
the Hubble parameter of the 4D space-time as a function of the density and to describe 
realistic cosmologies. It's interesting to study both analytically and numerically the system 
of Einstein equations. The analogous critical point analysis in the 5D bulk was performed 
in [43]. Some explicit solutions, derived with simplifying assumptions on the parameter of 
the internal space and on its geometry are also illustrated in our work. 

Further investigations on the 5D RS cosmology with brane-bulk energy exchange have 
been made from the holographic point of view [10], studying the dual theory in one lower 
dimension. The gauge/gravity duality [58, 59] (see [60] for a complete review) has undergone 
great improvements over the last ten years and provides a new approach to the analysis of 
brane-world models. As it is explained in [61]-[63] , the truncation of the AdSd+i space is 
equivalent to introducing in the dual picture a UV cutoff for the d-dimensional gauge theory 
(conformal field theory) coupled to d-dimcnsional gravity. Earlier suggestions about this 
idea are present in [64] . The presence of brane-bulk exchange corresponds to interactions 
between the gauge theory and the matter fields while the bulk "self-interaction" is shown to 
be related to the perturbation of the CFT (that becomes a strongly coupled gauge theory) . 
In [10] explicit examples of cosmological evolutions in the holographic 5D/4D picture as 
well as comparison between the two dual theories have been discussed. Other cosmological 
models have been analyzed in the context of the holographic correspondence [65]-[71]. 

Exploiting the AdS/CFT results, we build the holographic theory corresponding to the 
7D RS background. The 7D RS dual theory (see [10] for the analogous 5D/4D derivation) is 
then a renormahzed 6D CFT (the theory corresponding to the M5 system is an anomalous 
[72]-[76] (0,2) SCFT, but any other six dimensional large-N CFT can be chosen) coupled to 
6D gravity. See also [77] for other examples of holographic Weyl anomaly derivation. There 
are in addition higher order corrections to gravity and the six dimensional matter action. 
Higher derivative terms driven by conformal four dimensional anomaly [78] were proved to 
lead to an inflationary critical point in the 4D Starobinsky model [79] and to a successive 
greaceful exit from the long primordial inflation. As illustrated in [80] higher derivative 
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contributions to the Einstein equations cause the universe to enter a matter dominated era 
where the scale factor oscillates after inflation and to proceed through thermalization to a 
radiation dominated era. In our 6D holographic cosmological model we specially look for 
de Sitter fixed point solutions to the equations of motion describing late time acceleration 
or critical points suitable for early time inflation studying the associated stability matrix. 
A comparison with the 7D bulk analysis results shows some peculiar features of the 7D / 6D 
set-up. 

Interesting results for cosmologies with compactification emerge in the context of dynam- 
ical compactification [81, 82]. The compact space is treated with a different scale factor (as 
in the approach that will be used in this paper). In particular, in the context of dynamical 
compactification, the scale factor for the internal space has an inverse power dependence on 
the scale factor for the visible directions. The extra dimensions thus contract as the extended 
space expands. We also include some remarks on dynamical compactification applications 
in our set-up. Other attempts to reduce to conventional cosmology and investigate issues 
such as the cosmological constant from models with arbitrary number of extra dimensions 
are given in [83, 84]. In particular, cosmologies in six dimesions are analyzed in [85]. 

The structure of the paper is as follows. Next section will describe the set-up of the seven 
dimensional RS model. Section III will be focused on the cosmological evolution from the 
7D point of view, admitting brane-bulk energy exchange and particularly investigating the 
form of Einstein equations with some specific ansatz, while in section IV the critical point 
analysis will be illustrated, including numerical phase space portrait and explicit solutions. 
In section V and the following we will derive the 6D holographic dual to the 7D RS and the 
associated equations of motion. Section VII summarizes the fixed points in the holographic 
description and their stability. Some examples of the correspondence between the brane and 
bulk points of view will be given in section VIII. The generalization to non conformal and 
interacting theory, corresponding to non vanishing brane-bulk energy exchange and bulk self 
interaction in the 7D approach, will be exposed in section IX. Finally, the last section will 
summarize the results, give some conclusions and further considerations. The first appendix 
gives the form of the general anomaly for a 6D CFT in a curved space and of the other trace 
terms in the 6D theory and appendix B is devoted to the critical point analysis in the six 
dimensional cosmology. 
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II. 7D RS SET-UP 



As we announced in the introduction, we will work in a seven dimensional bulk with a 5- 
brane located at the origin of the direction z tranverse to the brane itself and with a z — > — 2; 
Z2 identification. In analogy to the 5D RS model, the action in this seven dimensional set- 
up is given by the sum of the Einstein-Hilbert action with 7D cosmological constant plus a 
contribution localised on the brane that represents the brane tension. Besides, we also put 
a matter term both in the bulk and on the brane. In formula we thus have 

•S" = SeH + SqH + Sm,B + Stens + Sm = 

= j d'^xV^ {M^R - A7 + + J d^xy/^ {-V + + Sgh (1) 

where V is the brane tension and we will call the associated contribution to the action 
Stens- Seh is the usual Einstein-Hilbert action with the seven dimensional bulk cosmological 
constant Ay and are, respectively, the bulk and brane matter lagrangians. Sgh — 

J d^x^J—^K, where K denotes the trace of the extrinsic curvature on the boundary, is the 
Gibbons-Hawking action added in order to cancel the boundary term arising computing the 
variation on the Einstein-Hilbert action and to get the usual Einstein equations. The action 
for the matter in the bulk Sm,B is an additional term with respect to the usual RS set-up, 
whereas the matter on the brane contribution will be refered to as Sm- The metric 7^,^ is 
the induced metric on the brane. The brane tension is necessary in the RS models in order 
to compensate for the presence of the cosmological constant in the bulk. 

The classical solution of the equations of motion for the theory above, neglecting all the 
matter terms and with a warped geometry of the kind ds^ = e^^dx^ + dz"^ {W = W{z) is 
the warp factor and the 6D x-directions are fiat), is the analogue of the solution described 
by Randall and Sundrum [1] for the 5D RSII model. The 7D solution gives as a result 
W{z) = 2| 

^ly ~ 30MS ' that the space— time is a slice of AdSj with the Z2 typical reflection, 
where there should exist a relation between the brane tension and the bulk cosmological 
constant 3V^ = -AOM^Aj. 

The aim of this section and of the next one is to generalize the RS ansatz to a time 
dependent background and to wrap the 5-brane over a two dimensional internal space, 
ending up with an effective 4D cosmology. Taking account of the seventh warped extra 
dimension and of the compactification over the other two extra dimensions, giving two 
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different warp factors to the 3D space and the internal 2D space, the time dependent ansatz 
for the metric is of the form 

ds^ = -n^(t, z)df + a^{t, z)Ciidx'dx^ + h^{t, z)iab^y''dy^ + f{t, z)dz^ (2) 

with the maximally symmetric Qj background in three spatial dimensions (with spatial 
curvature k) and $,ab for the 2D internal space (with spatial curvature k). We use capital 
indices A,B,... to run over the seven dimensions, i, j, ... for the three spatial dimensions 
of the 4D space-time, a,b, . . . for the two internal dimensions. In our notations z represents 
the seventh warped extra direction, the y coordinates belong to the 2D internal space, while 
the {Xfj_} = {t,Xi,ya} run over the 6D space-time on the brane. Summarizing, the structure 
of the bulk is thus made of a time coordinate, three extended maximally symmetric spatial 
dimensions (that gives, toghether with the time, the visible 4D space-time), two compact 
dimensions and a warped direction. The 3D and 2D spaces have two different scale factors 
a{t, z) and 6(t, z) respectively, while a gauge choice can be made for the values of the n(t, z) 
and /(t, z) factors on the brane, i.e. when z = 

A less physically meaningful background, but better understood, would be to have a five 
dimensional maximally symmetric space with some 5D metric and just one scale factor 
a{t,z), without compactifying on any two dimensional internal space. The solution to the 
equations of motion in this case is much simpler. The results related to this background will 
be briefly mentioned along with the more realistic analysis with brane wraping over the two 
dimensional internal space. 

III. COSMOLOGICAL EVOLUTION IN THE BULK 

In this section we will analyze some aspects of the cosmological evolution from the 7D 
bulk point of view. We will write the equations of motion for the bulk action and solve 
them by making assumptions to simplify their form and get explicit results evaluated on the 
brane. 

Given the set-up described in the previous section, we parametrize all the contributions 
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to the stress-energy tensor as 



TcU = ^ diag {-V, -V, -V, -V, -V, -V, 0) 



T^\v,B = diag (-A7, -A7, -A7, -Ay, -Ay, -A7, -A7) 
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Tc\m,b = -^diag(-p,p,p,p,7r,7r,0) 

Tc\m,B — Tq (3) 

with the subindices v and m indicating the vacuum and matter stress-energy tensors, while 
h and B stand for the brane and bulk contributions respectively. A difference between this 
(4+2+ 1)D background and the simpler (6+l)D analysis without the 2D compactification 
cited at the end of the previous section, is having in (3) two different pressures in the 3D 
space and in the 2D compact dimensions for the matter on the brane, while for the (6+l)D 
background we would put tt = p. This generalization is due to the fact that we don't assume 
homogeneity for the matter fluid in the whole (3+2)-dimensional space, but only in the 3D 
and 2D spaces separetely. 

Having calculated the Einstein tensor, we can put the explicit expression in the equation 

'^^^ = 2i^^^^ 

evaluated on the brane (from now on all the functions are evaluated on the brane, i.e. at 
z —>■ 0), m the specific background (2). As a consequence, for the 00, ij and ab components 
we obtain the jump equations 

K = = + ^ - 3P + 3^) (5) 

^ ^ 20M^ ^"^ + 4p + 3p + 27r) 

These are the values of the warp factors in the limit 2; ^ 0, where the subscripts + and 
— dinstinguish the limit taken from below from the limit taken from above. The prime 
denotes the partial derivative with respect to the z coordinate, while the dot indicates the 
time derivative. For the 07 and 77 components, substituing the expressions (5) and choosing 
a gauge with /(i, 0) = 1 and n{t, 0) = 1, we get the (non) conservation equation 

p + 3^(p + p) + 2^(p + 7r) = 2roy (6) 



and the Priedmann equation 



d b a? b"^ db k k 
3- + 2- + 3- + - + 6-- + 3- + - = 
a a'^ 0'^ ab b^ 

(6p - TT - p) + p (6p - TT + p) + + 



5^2 



+ ISA., - ^T/ (7) 



(20M 

+7 (P - tt) (y - 19p - 3p - Ttt) 

o 

We have defined the constant 

= 30^ ("^^ + ioIf^^O 
which plays the role of an effective cosmological constant on the brane. These (6)-(7) are 
two equations in five variables if, F, p,p, tt. We will thus have to make an ansatz for some 
of those variables. 

The pure RS system correspond to setting Tj = Tj = 0, that means putting to zero the 
brane-bulk energy exchange and no cosmological term on the brane, i.e. Xrs = 0, to restore 
RS fine-tuning. 

We can now write a simplified version of the differential equations (6) -(7), using the usual 
ansatz for the equation of state of the matter fluid on the brane, i.e. 

p^wp, TT = w^p (9) 



The set of equations, in terms of the Hubble parameters if = ^, F = |, is 



3iir + 2i^ + 6H' + 6HF + 3F' + 3^ + ^ = -j^ (cyV + c,p) p + - ^ (1°) 



p+[3(l + w)ii + 2(l + w^)F]p = 2To7 (11) 

with 



31^ — Gwtt — 5 llw + 14u'7r + 10 — — Wt^MZw — 7Wt^) 
= 400 ' = 400 



Looking at the definition of the two coefficients cy, Cp, we can note that equation (10) gets 
further simplified when the two pressures p and tt arc equal. This can be seen also from the 
previous equation (7), where the "non-standard" term on the r.h.s. (standard with respect 
to the homogeneous background analysis) is proportional to (p — tt) . We will first examine 
some cosmological solutions assuming p — tt and then we will drop this equal pressure 
condition to find an expression for H, in terms of the energy density, in the particular limits 
of static compact extra dimensions and equal scale factors. 



A. Equal pressures in 3D and 2D compact space 



We can first try to find an interesting solution by simplifying the computation assuming 
TT — p. In this case, the equation (7) written in terms of the Hubble parameters of the 3D 
space and 2D extra dimensions, defined as H — a/ a, F — h/h respectively, toghether with 
the (non) conservation equation (6), becomes 

[V^(5p-p)+p(5p + 2p)] + 15AK5-7^ (13) 



p+{3H + 2F){p + p) =2To7 (14) 

We note that the system of equations written above still contains three variables 
H{t), F{t), pit) but only two equations. So we are able to just determine the value of the 3D 
Hubble parameter H{t) as a function of the 2D one F{t). Moreover, given the complicated 
form of this set of equations, we will make some assumptions on the internal space, such as 
flat compact extra dimensions {k — 0) or static extra dimensions {F{t) = 0) in the following 
subsections. 

Manipulating (13) the system takes the form 

5^ {3H + 2Ff + 6 (3// + 2Ff + 6 (3i/ + 2F) {H - Ff = 
^ [Wp + 6 (3ii" + 2F) Vp + 5pp + 3 {2,H + 2F) p^] + 150 {3H + 2F) Xrs+ 

+ ^ iy + P) 27^7° - ]^ (3^ + 2^) - 10 (3iJ + 2F) ^3 A + (15) 



It is interesting to derive the first order ODEs from the second order one (15), in order to 
find the expression for H"^ as a function of the localized matter energy density p and to 
perform the critical point analysis. 



1. Flat compact extra dimensions with equal pressures 

In flat compact extra dimensions {k — 0) and flat 3D space (A; = 0), equation (15) shows 
that in the limit in which the two Hubble parameters are almost equal we can deduce the 
solution for (3if + 2F) in terms of the localized energy density p and of a mirage density x 
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that we will define below through a differential equation. In fact, in this case the third term 
on the l.h.s. of (15) is negligeable, leaving an analogous differential expression on both sides 
of the equality. 

The solution for {H — F) ^ {3H + 2F) and k = k, = is given by 

The solution is written in terms of the mirage density x the localized energy density p. 
The mirage density must satisfy 

X + ^ (3// + 2F) X = + l) 2r° -^(3H + 2F) T/ (17) 

and the equation for p is the (non) conservation equation 

p+{3H + 2F){p + p)^ 2To7 ( 18) 

Here wc get a linear and quadratic p dependence for the Hubble parameter H^, as well 
as a depedence on the mirage density x and on the hidden sector Hubble parameter F. The 
quadratic and linear terms in p are analogous to those in the 5D analysis [10], implying that 
for p ^ y the cosmological evolution looks four dimensional, while it moves away from the 
4D behavior for p^V. The term in x also already appears in the 5D model, as well as the 
Xrs constant term. Besides, the mirage energy density dynamics are controlled by the bulk 
parameters and Tj, that represent the brane-bulk energy exchange and bulk pressure, 
as in [10] . However, a new variable F, the internal dimension Hubble parameter, arises and 
remains undetermined unless we make an ansatz for it. We can argue that the solution (16)- 
(17) is written in terms of a "total" Hubble parameter ^{SH + 2F), that carries the same 
characteristics as the H Hubble parameter in the 5D model, but also includes the dynamics 
of the evolution of the extra dimensions. For equal scale factors, F = H, this "total" Hubble 
parameter reduces to H alone, giving the exact analogue to the 5D RS cosmology. 
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2. Equal scale factors with equal pressures 



A special case in which the {H — F) <^ {2>H + 2F) hmit is vahd is the equal scale factor 
case F — H. The results can directly be obtained from the previous subsection, yielding 

X + 6Hx=(f + l)2T»-l^i/r/ (20) 
/j+5ff(p + p) = 2r„7 (21) 

We added the curvature contributions that can be computed exactly in this limit. This 
solution is particularly simple thanks to the simultaneous vanishing of the {H — F) and 
(p — tt) terms. It shows the quadratic depandence of on p and the linear term in (p + x) . 
The mirage density reduces to free radiation in 6D space-time when we restrict to pure 
RS configuration, with no energy exchange. In this same limit, the localized matter energy 
density obeys to standard conservation equation in 6D. 

We will now drop the equal pressure ansatz and derive the expression for the Hubble 
parameter of the 3D space making particular assumptions on the internal space scale factor. 
We suppose from now on to live in a spatially flat universe {k — 0), where nevertheless the 
extra dimensions may be curved 7^ generally). 



B. Equal scale factors (generic pressures) 

For generic pressures, we make use of the parametrization by means of w,Wt^ for the 
pressures of both the non compact and internal dimensions p, tt that we introduced in (9), 
coming from the generalization of the equation of state for a fluid with energy density p. 

We evaluate the Priedmann and (non) conservation differential equations (lO)-(ll) as- 
suming the scale factors of the 3D space and of the 2D internal space to be equal and 
consequently assuming the Hubble parameters governing the cosmological evolution of the 
two spaces to be equal, F = H. The system (lO)-(ll) takes the form 

I + m') = {cyV + c,p) Hp-^ + IhHXns - (22) 

p+{'i{l+w) + 2{l + w^))Hp = 2TQ-r (23) 

with the coefficients cy, Cp still given by (12). With the help of (23), (22) can be brought in 
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a form from which we can exphcitly deduce H as a function of p and x 



+ 



yelding 



2T° T"^ k 

+ Mh (S^(-)'^ + 5p(e.)P) - + 15HXrs - -^H (24) 



^ - FTTTnP + r.^in (P + X) + Ars - — (25) 



5Mio^ 5M10 ■ ■ 10 a' 



X + = 2r° f 1 + 1^^) - :r^HT^ (26) 
p+{3{l + w) + 2{l + w^))Hp^2To7 (27) 



The two constants cvieq) , Cp(eq) must satisfy 

^ 3{l + w) + 2{l + w^)-6' ^ 3{l + w) + 2{l + w^)-3 ^^^^ 

in order to have the right coefficients in equation (24). For some values oiw,WT^ the denom- 
inator of cv(eq) or Cp(eq) may vanish. However we can fix such that cy (or Cp) becomes 
proportional to 3(1 + w) + 2(1 + w-,,) — 6 (or 3(1 + w) + 2(1 + Wt,) — 3 for Cp) and cv(eq) 
(or Cp(e,)) is finite. As an example consider = w (equal pressure in the internal space 
and 3D space, tt = p) and check that both Cv(eq) and Cp^^^^ remains finite and equal to 1/80. 
Clearly, when cv(eq) (or Cp(^eq)) diverges we cannot write the Friedmann equation (22) in the 
form (25). 



C. Static compact extra dimensions (generic pressures) 

We can follow the same procedure as in the equal scale factor hmit for the case of 
static compact extra dimensions F = 0. While in the previous subsection the two internal 
and observed spaces were evolving according to the same dynamics, in this limit the extra 
dimensions do not evolve and remain static. 

The two differential equations of motion for the gravity action in this case become 

p+3{w+l)Hp^2To7 (30) 
where cy, Cp are as before (12). We introduce the new coefficients cv{st)-i Cp(^t) defined by 
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After plugging (30) into the Priedmann equation (29) we come to the expressions for H and 
X 



X + 4i^x= (l + ^^)2T°--^i^r/ (33) 
p + 3{w + l)Hp^2Tor (34) 

In analogy to the equal scale factor case, these expressions are valid as long as we don't have 
w = 1/3 {w = -1/3) with cvy^O {cp ^ 0). 



D. Proportional Hubble parameters 

We are going to combine in the same description the two hmits of a(t) = 6(t) and -F = 0, 
impliying, in the first case, an equal cosmological evolution for the internal space and the 3D 
visible spatial dimensions and, in the second case, the absence of evolution for the compact 
space. 

Both the two systems of differential equations (25)-(27) and (32)-(34) obtained in the two 
different limits can be written in a unified formulation that encloses the two just cited sets 
of equations, defining some appropriate constant parameters. We introduce an "effective" 
number of dimensions d that takes the values d = 6 in the equal scale factor limit and = 4 
in the static compact extra dimensions. If we look at the equations (26) and (33), we see 
that d appears as the number of dimensions for which the energy density x satisfies the 
free radiation conservation equation in the limit of pure RS (Ty = Ty = 0). In fact, the 
Priedmann equation plus the two (non) conservation equations can be rewritten as 

^ ^(rf-i)Mio^ ^ {d-m^^^p^^^'W^)'^^ W^) ^ ^ 



p + WdHp^2To'j (37) 

We have in addition defined = 3(1 + w) + (c? — 4)(1 + Wt^) and cv,d — c-vlk^d ~ d), Cp^d — 
Cp/ {wd — d/'^) where cy and Cp are given in (12). We remind that in order to get an algebraic 
equation for as a function of the energy densities p and x (35), we have to keep cv,d-i Cp^d 
finite, i.e. respectively Wd ^ d,Wdy^ d/2 unless cy = 0, Cp = 0. For example we cannot write 
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H in the form (35) ii w — l/3,w.,r — in both the equal scale factor and the static compact 
extra dimension limit, since the linear term in p has a diverging coefficient. 

Moreover, we can futher generalize this analysis introducing a parameter ^ such that 
F = ^H. This description contains all the above studied limits. The analogous of previous 
relations (35)-(37) can be written as 

(3 + 2e)Mio^ (3 + 2e)Mio + 3^ + 6 a?^ + 2^ + 2 

(38) 

X + ^./^X = 2To(l + |i^)-^Hr/ (39) 
p + w^Hp = 2ro7 (40) 

where now is a more complicated function of the proportionality constant ^ between the 
two Hubble parameters = 6^^^^^ and it reduces to d = 6, d = 4 in the two previously 
examined limits of equal scale factors and static compact extra dimensions (^ = 1,^ = 0). 
The constant reduces to for { = 0, = 1 and is defined by W(_ = 3(1 + w) + 2^(1 + w^). 
The two coefficients Cy,^,Cp^^ arc defined with Wd — > wg, d d^. The result (38) 

is valid unless ^ — —3/2. In that case the equation for H becomes algebraic — though we 
still have the curvature term explicitely depending on the scale factor — , thus 

rr2 i 2 ^ a \ /A-t\ 

^ = -SSM^oP - - ^ + 4^^^ (41) 



p + 3{w — Wt^)Hp — 2Tc 
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No mirage density appears and the Hubble parameter is a quadratic polynomial in the 
localized energy density p alone. Besides, if the pressures are equal Wj^ = w in the pure RS 
set-up Ty = 0, the energy density is constant in time and so is H'^ + for ^ = —3/2 (41). 
The set of equations (41) also doesn't depend on T/ at all. 

Again we have to restrict to ^ d^,d^/2 to keep cv,^, Cp^^ finite (unless cy oc w^ — d^,Cp oc 
w^-d^/2). 

We remark that for the scale factors satisfying b{t) = l/a{t), i.e. dynamical compacti- 
fication with ,^ = — 1, the equation for the mirage energy density x in the pure RS set-up 
is still an effective 6D free radiation conservation equation, as for b{t) = a{t). In fact, the 
only solutions to = 6 are ^ = ±1. To obtain a 4D free radiation equation for x we have 
to require ^ = 0, since the second solution to = 4 is ^ = —3/2, for wich we don't define a 
mirage density (41). 
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However, this is not the end of the story. Introducing the effective 4D densties g — V(2)P 
and X = V{2)X (where V(2) = vb'^{t) — va'^^{t) is the volume of the 2D internal space), 
we have to replace the l.h.s. of equations (39)-(40) repsectively by ^ + {d^ — 2{)i?;t and 
g + {w^ — 2^)Hq (the r.h.s. are also modified and we will explicitly write them at the end of 
subsection IV B). This tells us that the 4D mirage density is a free radiation energy density 
for pure RS in four dimensions for — — 4, which has solutions ^ = 0, ^ = 1 — i.e. static 
internal space or equal scale factors, justifying the study of these two limits. 

E. Comments 

We here summarize some considerations about the bulk evolution equations derived in 
the previous subsections and, in particular, about the explicit expressions we have found for 
the Hubble parameters in the discussed limits. 

- With the assumption of having the same pressure for the matter fluid in the two 
dimensional internal space and in the 3D visible space (i.e. tt = p), we found a form of 
the Friedmann equation that has the advantage of keeping both the Hubble parameters 
not constrained by any particular ansatz. The Friedmann equation (15) provides an 
expression for (3i? + 2F) in terms of p and x (fo^" spatially flat spaces). This solution, 
though, is satisfactory only in the hmit of small {H — F). When [H — F) is not 
negligeable w.r.t. (3i7 + 2F), the mirage density equation may be written introducing 
an extra term independent of the bulk parameters Tj,Tj. This prevents x to obey to 
a free radiation equation in the pure RS set-up (T^ = Ty = 0), as it should instead 
be in the context of the AdS/CFT correspondence (we will discuss the comparison 
between the bulk and the dual brane analysis in section VIII). 

- In the simple limits of equal scale factors (25) and static compact extra dimensions (32) 
we recovered an expression for containing a quadratic term in p and a linear term 
in (p + x) ) where p is the localized energy density and x is an artiflcially introduced 
mirage density that accounts for the bulk dynamics. In fact it depends on the bulk 
parameters Ty , Ty . When Ty = = 0, the mirage density obeys to 4D free radiation 
equation for the static compact extra dimension case and to 6D free radiation equation 
for the equal scale factor case. This is in complete analogy to the 5D analysis [43], 



16 



where the same dependence on p and x occurs and the mirage energy satisfies 4D free 
radiation for pure RS (i.e. = r| = 0). 

- When — w (or equivalently tt — p) in the equal scale factor limit (25), we find the 
results given by the equal pressures subsection in the case F = if (19). The two hmits 
of equal pressures and equal scale factors then commute and the results are consistent. 

- The description of section HID encloses in a unifying way the results in the limits 
of static compact extra dimensions and equal scale factors. It moreover generalizes 
these results to the case of evolutions governed by proportional Hubble parameters 
F{t) = ^H{t). We will use the set of equations written in terms of the effective 
number of dimensions d (that describes the two limits of static internal space, with 
d — 4, and equal scaling for the compactification space and the 3D space, with d = 6) 
to study the corresponding cosmological evolution in the next section. 

We are now going to proceed to the analysis of the critical points for this seven dimensional 
universe in a 7D Randall-Sundrum set-up, including the energy exchange term. 



IV. BULK CRITICAL POINT ANALYSIS WITH ENERGY EXCHANGE 



We have until now transformed the second order differential equation (7) plus the (non) 
conservation equation (6) in a set of three linear differential equations (35)-(37) for combined 
equal scale factor and static compact extra dimension limits, or more generally (38)-(40) 
for proportional Hubble parameters. We have introduced the mirage density x defined by 
its differential equation. In this section we will use the first system of equations (35)- 
(37), obtained to describe both the hmit of equal scale factors and static compact extra 
dimensions, to find its fixed points and the corresponding stability. The critical point analysis 
will allow us to study the cosmological evolution in the bulk description for F = or F = if. 

We make an assumption on the bulk components of the stress-energy tensor that appears 
in the differential equations for the energy densities p and x- As in [43], we will take the 
diagonal elements of the stress energy tensor to satisfy the relation 



(diag) 



m,B 
p{diag) 
-v,B 



(diag) 



m,b 
p(diag) 
" v,b 



(42) 
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This enforce the solution to the Priedmann equation to be reasonably independent of the 
bulk dynamics, since the Ty term in (7) becomes negligeable with respect to the first term 
on the r.h.s. of the same equation. Imposing such a relation, T/ disappears from the sets 
of linear differential equations, while we remain left with the component. For the future 
bulk calculations we will define T = 2T° to simplify the notation. 

Before starting the critical point analysis we note that when T — the system of equations 
(35)-(37) have only trivial critical points characterized by zero visible Hubble parameter 
when the internal space is flat. There are two of these critical points. One is given by 
= —K/2{d— 1)6^, = X* = (which is vahd only if we are compactifying on hyperbolic 
or fiat spaces) and the other is = 0, Cp^^pl + 2cv,d^ {p* + X*) = M^^n/h. 

We will first restrict to small density [92]p ^ V and flat internal space k = (remind 
that the 3D space is already supposed to be flat, having put /c = 0) and then go through the 
generic density analysis. Eventual de Sitter stable solutions (for the 4D visible space-time) 
could represent the present accelerated era, while inflationary phases at early times may be 
associated to primordial inflation. 



A. Small energy density and flat compact extra dimensions 

When the localized energy density is small and the internal space curvature vanishes, 
p <^ y, K = 0, the bulk Einstein equations (35)-(37) in terms of x P become 

^'-i^oiP^x) (43) 
X + dHx^T (44) 
p + waHp = -T (45) 

We note that in this approximation expression (43) isn't valid for — d, unless cy oc {w^—d) 
(this would for example determine a specific value for w^^ as a function oiw). Nevertheless 
we will find in the rest of the section that cy.d always appears with the coefficient [wd — d) 
in the critical point analysis. We thus have to keep in mind that (w^ — d)cv,d = Cy ~ 
(31w — 6w,r — 5)/400 always remains finite. 

Fixed point solutions In the small density approximation, the fixed points in terms of 
the critical energy density can immediately be found (to have a full solution we have to 
make an ansatz on the form of the energy exchange parameter) . 
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The first solution is given by 




(46) 



d{d- 



<)'2cv,dV \ 
l)Mio J 



1/2 



where = —wa 



This represents an inflationary critical point for the 



cosmological evolution, for {wd — d)cv,d — cy > 0, i.e. both cv,d > and d > Wd or cv,d < 
and d < Wd [93]. The acceleration factor at the flxed point is simply given by q^, — H^. 
Since we assume Wd to be positive, is negative and H^, in (46) is positive. At this flxed 
point the universe is thus expanding. 

Another fixed point leaves x* unchanged, while and have switched signs with 
respect to (46), meaning that is negative and the universe in contracting (there is a 
symmetry T ^-T,H ^ -H) 



The trivial critical point is characterized by mirage density equal and opposite to p^, but 
zero Hubble parameter and energy exchange (in the case we admit for the energy exchange 
the form T = Ap" all the variables are zero at the trivial fixed point). 

For positive critical energy densities p^, we obtain a negative brane-bulk energy exchange 
parameter at the critical point if if^ > and, viceversa, wc have positive critical energy 
exchange for a contracting universe at the critical point. During the evolution, we can 
expect a change of regime going from negative to positive T as the energy density localized 
on the brane grows, as for the 5D RS critical point analysis with energy exchange carried 
in [43]. Even though most of the analysis in this section will be performed supposing that 
the energy exchange parameter has fixed sign determined by the sign of A (since we will 
mainly assume T = Ap'^), we can argue that for small energy density p the generic energy 
exchange is presumably negative, meaning that energy would be transferred from the bulk 
onto the brane. In this hypothesis, an equilibrium can be reached, such that the energy 
density would have a large limiting value for which energy starts to fiow back into the bulk 
(with positive energy exchange). 




(47) 
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Stability analysis The real parts of the eigenvalues of the stability matrix for the {5x, 5p) 
linear perturbations corresponding to the critical point (46) can have opposite signs or be 
both negative. This depends on the value of T as a function of the energy density p at the 
fixed point (46) describing an expanding universe with energy influx. The explicit form of 
the eigenvalues is 

A± = ^pI'^ \d+{l- v)wd ± ^{d + (1 - v)wdY - 2(3 - '2D)dw^ (48) 

ZWd L -I 

where we have defined 

aiog|T 



V = 



(49) 



9 log p 

Expression (48) then shows that the two eigenvalues have negative real part when z> < 
3/2. There is a second upper bound on i> derived from requiring negative real part for the 
eigenvalues. Nonetheless, for the range of values — 1 < tf, < 1 and d = 4, 6 in which we're 
interested, this bound is always equal or greater than 3/2. If \T\ is a decreasing or constant 
function of p near T^,, the non trivial inflationary critical point always is an attractor. Also for 
growing |T| we can have stable inflationary fixed points, as long as z/ < 3/2. In particular, 
the linear u = 1 case is included in the stable inflationary fixed point window and will 
be analyzed both solving the Einstein equations numerically, in the next subsection, and 
deriving an explicit solution, in subsection IV C. 
Besides, when 

/ 2d d _ ^ f2d d . 
1-W <i/<l + W (50) 

\ Wd Wd V Wd Wd 

the eigenvalues have non zero immaginary part and the critical point will hence be a stable 
spiral for 9 < 3/2. For values of v out of the range (50), we get a node. As an example, 
let's assume the value v — I'm the equal scale factor background. This gives a stable spiral 
for > 1/2 or < 1/2 and w > —2(1 + Wtt) (considering w, > —1). This means that 
in the case -u; ^ 1/3 and w^^ = the critical point is a stable spiral, while for both w and 
null we instead have a stable node. 

For energy outflow > (which goes along with contraction < 0), we get a minus 
sign overall modifying the eigenvalues (48) referring to the linearized system around the 
(47) critical point (characterized indeed by energy outflow). The eigenvalues cannot be 
both negative in this case. In fact we should demand Wd > d/{i> — 1) with u > 1 but also 
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u < 3/2 to get a stable fixed point. Only the trivial point, as we will discuss later, can be 
attractive for energy outflow dynamics. 

Assumption T = Ap^ and numerical solutions Assuming the brane-bulk energy ex- 
change parameter to take the form T = Ap'^ (so that v = v referring to (49)), we can rewrite 
the system of differential equations in term of dimensionless quantities p = 7^p, x = l^'^j 
H — ^H, T — 7^T, where we called 7^ = -pzlpjiu-- The dimensionless variable p/V used to 
perform the small energy density expansion at the beginnning of this section is related to 
the dimensionless variable phy p — (^cy^J{d — 1))^^^ (y/M^)^^^ {pI^)- So, considering the 
small p/V approximation is equivalent to considering small p approximation if the brane 
tension V satisfies V < with respect to the 7D Planck mass and cy^ is reasonably of 
the order cv,d ^ 1- The complete set of fixed point solutions (discarding the trivial ones) 
can be calculated in terms of the parameters A, v characterizing the energy exchange, wa, d 
denoting the background (static extra dimensions or equal scale factors) and the equations 
of state for both the 3D and internal spaces. 

The Einstein equations become 

H'^ = cv,d (p + x) 
jc + dHx = Ap^ (51) 
p + WdHp = -Ap"" 

where A = 7^+6(^-^)^. 

The acceleration q can be evaluated independently of 1/ and A 

^=(^~t)^^"^'^+O~0^^"'^ ^^^^ 

as a function of the localized matter density and of the mirage density. Due to the positive- 
ness constraint on cy^^ (p + x) coming from the first equation in (51), the trajectories in the 
phase space must satisfy 

q < {d - Wd)cv,d P (53) 

as it is indeed showed in the numerical plots of figure 1. For > 2 and cv,d > we 
have positive acceleration only if the mirage density x is negative and smaller than — p(2 — 
Wd)/{2 — d). On the other hand, x gets positive (suppose cv,d > 0) only if g < (2 — Wd)p/2. 
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The fixed points are given by 



X. 



Z-2u 



\3-2u 



\3-2u 



d 

^2(.-l)^2 



Wd 



P 



~ 3-2u 



(54) 
(55) 
(56) 



wlcy^did - Wd) 

For V < 3/2, when the non trivial fixed point is stable, we have two roots of (56) with 
opposite signs if z/ = 1/2 + m,m G Z. Only one real root exists for integer v and it carries 
the sign of the r.h.s. in (56). Finally, for v = (2m + l)/4, we have a positive root if the r.h.s. 
in (56) is negative. The two eigenvalues corresponding to a negative p^, always have opposite 
real part, implying that this fixed point isn't be stable at linear order, it is a saddle. We 
further note that for Wd > d and positive cv,d (or alternatively d > Wd and negative cv,d) we 
can only have real and positive fixed point if z/ = (2m + l)/4, m G Z. This implies that 
for u = 1 there is no non trivial fixed point with positive p^, when {d — Wd) and cv,d have 
opposite signs. Moreover, these negative p^, points are always characterized by negative H^, 
so that they wouldn't be infiationary. 





(a) (b) 

FIG. 1: Phase spaces q/p for different initial conditions pQ and %, with A = —1, v = 1, d = 6 
(equal scale factors — analogous pictures come from the static extra dimension case) and: (a) 
Wd = 4: (for istance w = OjWtt = —1/2 or w = —1/5,11!^ = 0) leading to a spiral like stable critical 
point, (b) Wd = 3 (for istance w = OjWj^ = 1/2) determining a stable node. The extra blue grid 
lines intersection represents the fixed point. 



We can check the stability of the critical points by means of a numerical analysis of the 
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differential system of equations (51). In the case of energy influx A < 0, putting u — 1 and 
different values for the d,Wd parameters, we get the phase spaces in figure 1, plotting the 
acceleration factor q{t) = a(t)/d(t) as a function of the energy density p(t). We thus check 
that, solving the system of differential equations for variable initial conditions for x P-, 
including both positive and negative initial q, all the different trajectories converge to the 
non trivial fixed point, designated by the intersection of the two perpendicular lines in the 
picture. Besides, as we expect, in part 1(a) they have a spiral behavior, while in the 1(b) 
case they denote a node. 

In the limit Wd ^ dwe numerically recover the analytical solution discussed in subsection 
IV C, neglecting the very large density behavior. 

For u > 3/2, i.e. when the non trivial fixed point is no more an attractor, we find that 
some of the trajectories go to the trivial critical point, while another branch of solutions 
to the Eintein equations (51) are characterized by diverging energy density p (they become 
unreliable when > {2cv,d/{d— l))^{V/M^y). This happens because, as it is suggested by 
the integration of the third equation in (51) with energy influx hypothesis, for p big enough 

~ 1 /2 

— precisely for > (w^/A) Cy^ (1 + x/p) — the function p{t) starts growing, while 

for small p, p{t) eventually goes to zero. Depending on the initial conditions we will have 
solutions ending in the trivial fixed point or diverging. 

The behavior of the system with energy outflow can be analytically deduced. Given 
the hypothesis T > it is clear that the trajectories can't be attracted by the fixed point 
solution — B^,pl^^. They may go to the critical point characterized by negative Hubble 

3 /2 

parameter and = —B^pJ . However, we already determined the non attractive nature 
of this fixed point. Another way to describe outfiow dynamics is to conclude from the form 
of the Einstein equations (45) that all the trajectories in the phase space q/ p go toward the 
trivial point, since for positive T the density p is suppressed at late time. The way in which 
the trajectories go to the critical point depends as before on the positiveness of the function 
under the square root in (48). We numerically checked that ior u — 1 and d — 6,Wd — 4: the 
null fixed point is an attractor and, in particular, a stable node. 
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B. Critical points with general energy density 



Allowing K to be different from zero, not restricting the localized energy density to be 
small and with the assumptions (42) on the bulk matter stress-energy tensor, we have to 
solve the following set of equations 

Tcv.V , . 1 



-"PA 



{d-l)M 



10 



(ci- l)Mio 



(p + x) 



K 



2((i- 1)62 



\ cv,4 V 
p + WdHp = -T 



(57) 

(58) 
(59) 



As in the small energy density regime, we can immediately notice that the analytical 
behavior for energy outflow is characterized by decreasing p in time. The trajectories will 
thus be attracted to the trivial fixed point characterized by vanishing p^. 

We will now carry a general critical point analysis. Equation (57) exibits divergences if 
Wd = d Wd = d/2, unless cy oc {wd — d) or Cp oc {wd — d/2). When the divergence arises, 
the right system of equations is (10). Terms like {wd — d)cv,d — cy and {wd — d/2)cp^d — Cp 
are always finite (where we recall cy = 31 w — Owtt — 5 and 0^ = 11^ + 14u'7r + 10 — (w — 
wt:){?)W — Iwt^)). 

Fixed point solutions If we demand H to be positive, i.e. expanding universe, the 
solution for and is given by 



X* ~ 



-^P'J' 



Wd 
Wd 

d 



1 + ^ 



Cp^dP-k 



(60) 



b.pT 



where = 5* (p*) depends also on k and and is defined by 

{d - 2w)cp^aP* + 2{d - w)cv,dV 



bAp^ 



-Wd 



1/2 



(61) 



d{d-l)M^o 2{d-l)blp^ 

We have a negative energy exchange parameter, as in the small density limit, and a positive 
Hubble parameter. 

The second fixed point solution is equal to the first except for the Hi, and signs reversed 
(keep in mind the H — > — i7, T —>■ —T symmetry), such that would be negative and we 
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would have energy outflow at the critical point 

r. = B,iST 

The trivial critical point is characterized by vanishing i?* and Ty^, while the mirage density 
becomes 

X* = 7I?^- l + TTF^ (63) 



If the energy exchange is supposed to be of the form T = Ap", the trivial fixed point 
is characterized by zero value for all the variables except for the mirage density x which 
becomes x* — M^'^K/4blcv,dV and is zero for flat compact spaces. 

As in the limit of small energy density considered in the previous section, the constant 
is negative whenever the argument of the square root in (61) is positive, i.e. when 

K {d - 2wd)cp^dp^ + 2{d - w)cv,dV 

< jn — IX ,^in — P* (64) 



2{d-l)bl d(d-l)Mio 

If the square root gives an immaginary number, we don't have any real valued fixed point 
except for the trivial one. 

Stability analysis The positiveness of the eigenvalues of the stability matrix depends 
now on all the parameters and constants of the theory and not only on i>, as for the small 
density, flat compact extra dimension simple case. We consider the situation in which the 
variation of vanishes (this happens in the static compact extra dimension limit or for 
K = in the equal scale factor background, as the ratio remains constant). Otherwise 
we would have a linearized system of two differential equations plus one algebraic equation 
in the four variables 5k, 5p, 5x, 5H. 

There are two conditions that must be satisfied, in order to get two negative eigenvalues 
and hence a stable fixed point. These conditions give two upper bounds for i> in terms of 
the constants Wd, d and V, p^, 

~ T (d - 1)!b2M^o ~ '^Wd)cp,dPi. + (d- Wd)cv,dV] (65) 
1/ - 1 < — 

Wd 
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The second bound in (65) satisfies d/wa + 1 > 3/2 in the range —1 <w,Wt^< 1. Besides, 
the first bound reduces to v < 3/2 when we take the hmit p/V and put k — Q. So the 
results are in agreement with the previous small density analysis. 

The bounds (65) depend on the fixed point value of p, which can't be determined without 
making any assumption on the form of T. However, we can make some remarks on the nature 
of the fixed points. For values of D in the range 

1- — - i?*<i><l- — + i?* (66) 

Wd Wd 

where we defined 

the stability matrix eigenvalues have non null immaginary part and the trajectories near to 
the critical point have a spiral hke behavior. When i?^ < we always have node like fixed 
points. In agreement with the small density case (50), when p/V -C 1 and k = we get 



R^^^/2d/^d■ 

As an example, we assign the value i> — 1. Since the first bound (65) can be rewritten as 

this means that we should have i?^ > to get stabihty. We also find that the fixed points 
have spiral shape when > d/wd or i?* < —d/wd, they will be nodes otherwise. 

Assumption T = Ap^ and numerical solutions To do a more quantitative analysis we 
have to make an ansatz on the form of the energy exchange parameter T. As in the previ- 
ous section, we suppose a power dependence on the energy density p such that T — Ap'^. 
The equations for generic energy densities and internal space curvature can be rewritten 
introducing dimensionless variables as in (51) 

H"^ = Cp,d ap^ + cv,d (p + x) - K 



X + dHx = Ap''{l + 2^ap (68) 
V Cv,d ) 

p + Wdiip = -Ap" 

where a is a dimensionless constant defined by o? = -^^^j^ (t") ' ^ dimensionless 

2 

variable k — 2{d-i)b^ — ^® remind that we restrict to constant k approximation. 
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To obtain real observables, we have to restrict the possible values for p and x such that 
Cp,d otp^ + cv,d (p + x) — > 0. In fact, the plots show the presence of a prohibited zone 
in the phase space — in particular, in figure 2(a) it is clear that the region of the possible 
trajectories is delimited by a parabola. The relation that must be satisfied, in terms of the 
acceleration parameter and the energy density, is 

q <{d- 2wd)Cp^d otp^ + {d- Wd)cv,d p - dn (69) 

In fact, the analytical expression for the acceleration q — H + can be written using (68) 
in terms of the visible energy density and the mirage density. For any v 

? = (1 - Wd)cp,d ap^ + (l ~ ■y ) ^^.-i ^ + ("^ ~ ^^^^ 
So, taking a specific value for p, we can have positive acceleration for our universe only if 

- {Cp,d + cy,d P - < cv,i X < ^^^) 

and a necessary condition for this to be possible is a bound on the energy density [d — 
2wd)cv,d(^p^ + (o? — Wd)cv,dP > dk, as we can deduce from (69). The mirage density x has 
to be negative to get positive acceleration for > 2, k > 0. If instead Wd < 1, k < 0, the 
mirage density is positive for negative q. 

Manipulating the set of equations (68), we write the following differential equations in 
terms of the generic energy-exchange parameter T 

= -dx + r]T (l + 2^ap] [cp,^ ap^ + cv,d (p + x) - k] (72) 
= -WdP - r]f [cp,d ap^ + cv,d (p + x) - ~' (73) 
We thus come to the differential equation for the acceleration factor 



+ 



2q;(1 - Wd)cpp+ i(2 - Wd)cv 



p + dq (75) 



where H = {2aCpp^ + cyp + dk) /2 + q and rj = ±.1 denotes the two possible roots for 
H. Again we note the presence of the symmetry H —H, T —T. We have used the 
definitions {wd — d)cv,d — cy and {wd — d/2)cp^d — Cp- Prom this equation we can infer that 
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positive q implies growing q in an expanding universe (77 = +1) with energy outflow (T > 0) 
if cy < 0, > and p < fj^i^ = Pum- For cv,d: Cp^d > 0, cy < 0, Cp > this is reahzed if 
d/2 < Wd < d. In the case of energy influx (T < 0), we get increasing positive acceleration 

if Cv,Cp > and < 1 for all positive energy densities (cy^diCp,d has to be negative). Or 
else, q grows as p grows if cy > 0, Cp < and wa < 1, until the energy density reaches the 
bound piim- 

The non trivial flxed points for energy influx are determined by the roots of the equation 
for p^ 

1 1 

- ((i - 2wd) Cp^dOipl + -zcv,d (d - Wd) p* 2 ^^^"""^^ - K = (76) 

while for and we get the two functions of 



d \ cv,d J Wd 
Wc thus have to fix a particular value for v in order to establish the precise number of roots 
and the explicit solution for the critical points. For integer and seminteger v the number of 
roots we can obtain, keeping d ^ 2wd, is 

u>2 =^ # roots = 2(1/ - 1) > 2 
1 < i/ < 2 =^ # roots = 2 

i/ < 1 =^ # roots = 2(2 - I/) > 2 

These are all the roots of (76), including trivial and complex roots. When d — 2wd the 
critical value for the mirage energy density diverges due to the divergence of Cp^^, unless we 
fix to keep it finite. In this case, the number of roots changes to 

u>l =^ # roots = 2(i/ - 1) > 1 

p = 1 =^ # roots = 1 

v <1 ^> # roots = 3 - 2i/ > 1 

For 1/ > 1, one of the roots of equation (76) is null if k — 0. 

There are moreover two trivial flxed point solutions given by p^ — — 0, = cv,d ^ 
and Pi, — x-k — 0, Hi, — \/—k, that reduce to a unique point with all vanishing variables 
when the internal space is flat. 
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Let's study in more detail the case i' — 1, since solutions can be written explicitely being 
the case with the minimum number of roots for (76), toghether with the v — 2 case. As 
a result we get a trivial fixed point solution with = = 0, x* = cy^d and the trivial 

solution, acceptable only for negative and zero curvature, H^, = \/—k, = = 0. Finally, 
the two non trivial solutions (for d ^ 2wd) are given by 



-cv,d{d - Wd) - 4(d - 2wd)cp^d OiK'^ ± ^ [d - Wd^c^^^ - ^^i^ - '^Wd)cp^d aK^ 
^ 2id-2wd)cp,,a ^^^^ 



-cvAd - Wd) ± {d - Wd)'^c^^^ - 4:d{d - 2wd)cp^d otK'^ 
^* 2{d-2wd)cp^aOi 

where K corresponds to a shift and rescaling of due to the nonvanishing value of k and is 
defined by K'^ = /w^ + k. The two roots are both real only if the argument of the square 
root in (78) is positive, i.e. when Wd lies outside the two roots —ad {l ± \J {a + l)/a^ , 
with a = Acp^^aK — 1. If a is in the range bounded by —1 and the square root is 
always real, whatever d,Wd we choose. If two or all among cv,d,Cp,a, {d — Wd), {d — 2wd) 
have equal sign, one of the two solutions (78) always is characterized by a negative p*. 
We note that we can have at list a non trivial fixed point with positive energy density 
if {d — 2wd)cp^dO( > 0,{d — Wd)cv,d < — exactly two positive p^ fixed points — , or for 
{d — 2wd)cp^d a < — only one critical point with positive energy density. 
The non trivial solution for d = 2wd is 

H^^- — , p^ = 2cy,,X^ x*^-cv,dK^Acp,,K^a + l) (79) 

Wd 

For this unique fixed point solution to be characterized by positive p^, we have to demand a 
positive cv,dK'^- Moreover, we can derive from (67) that for k > cv,dWdP*/4:d the fixed point 
is a spiral, so that for istance, in a flat internal space, we always obtian a node since cv,d 
must be positive in order to have a positive fixed point {K = A in this case). 

The numerical analysis can now show some of the features that we commented for the 
cosmological evolution with generic density. The differential system of equations (68) (sub- 
stituing some precise values for i/) can be solved numerically in order to check the existence 
of stable inflationary critical points. In figure 2 we plot the dimensionless acceleration factor 
q{t) as a function of the dimensionless energy density p{t), as we did in the previous section 

29 



for small densities. In the plots, we consider for simplicity positive cv,d-, Cp^d and flat compact 
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FIG. 2: Trajectories in the phase spaces q/p with varing po and XO) d = 6, A = —1 (analogous 
pictures come from the static extra dimension case) and: (a) = 4 > d/2, a = 1, z/ = 1 leading to 
stable node in = 1/4 plus a second repulsive node in = 3/4, (b) Wd = 2.5 < d/2, a = l^u = 1 
determining a stable node, (c) Wd = 2.5 < d/2, a = l,z/ = — 1 an example of a stable spiral, (d) 
Wd = 2 < d/2, a = 10^, V = 1 sl spiral behavior for large a (large /V). 



extra dimensions k = (k 7^ results in a shift for A^/w'\ in the critical point evaluation 
and in some scaling of the suitable value of v in order to have stability). 

All the critical points we get are characterized by positive q, i.e. they represent an 
inflationary point. In the phase space portrait 2(a) trajectories starting with positive accel- 
eration factor and energy density lower than the critical one pass through an era of larger 
acceleration and then slow down to the fixed point where we have inflation. There are then 
solutions starting with negative acceleration and going to the positive q critical point, even- 
tually passing through a larger acceleration phase or through a smaller density phase. The 
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families of solutions that distinguish the diagram 2(a) from the others are caracterized by 
both initial and final very high energy density, since they are repelled by the second non 
attractive fixed point. Some start with high energy density and negative acceleration at 
late time, go through an era of larger acceleration (eventually positive) and then, while p 
becomes very large, they go to a region of large and negative q. Other go from positive 
aceleration to large negative q and large p. The mirage density x can start from an initial 
condition smaller or bigger than the critical value, has to be positive for g < — (3p + 1) p 
— as we deduce from (70) plugging in the values for the parameters — and approaches a 
negative constant value. 

In diagram 2(b), trajectories starting with negative acceleration go to the positive q 
fixed point, eventually reaching a maximum q before ending into the critical point. Positive 
acceleration initial condition lead to growing acceleration at very early times, when usually 
the energy density grows as well, then both q and p decrase to reach the fixed point, passing 
through a minimum for acceleration. The mirage density goes to the negative critical value 
being initially positive for the trajectories that come from negative acceleration conditions, 
with q < — (4p + l)p/2. With the choice of parameters we used, we get Wd < d/2. We 
could also have used the fixed point solution (79) if Wd = d/2 = 3 and the diagram for the 
phase space would have been analogous to that in plot 2(c), keeping the same values of plot 
2(b) for the other parameters. 

In the phase space 2(c) all solutions converge to the fixed point with a spiral behavior. 
Energy density and acceleration parameter thus oscillate around the critical values. Here 
u = —1 and d < 2wd. The number of roots corresponding to the critical point solutions 
are six, but four of them are complex roots and one is characterized by negative energy 
density. Only one real critical point with positive energy density exists in this case and it 
has Hubble parameter. Since q can be negative at some time of the evolution, even if starting 
with a positive value, x can pass through a positive phase, crossing zero before reaching the 
negative critical value, if ^ < — (4p + 1) p/2. A similar plot can be drawn also if Wd — d/2. 
There would be only five roots for p^,, four of which would be complex conjugated and the 
last would have positive energy density, representing the stable spiral. 

Another stable spiral is represented in figure 2(d). Here, the dimensionless parameter a, 
which is proportional to M^/V, is large and u — 1. We find only one non trivial fixed point 
with positive energy density and spiral behavior, so that trajectories has a shape analogous 
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to the ones in 2(c). 

For values of v different from v — 1 the number of fixed point roots may vary according 
to the previous discussion. Nonetheless, (as it is shown as an example in figure 2(c) for 
V = —1) some of the roots may be complex conjugated and thus not acceptable. Another 
simple case is z/ = 2, where we get two solutions to (76), as with the v = 1 assumption. We 
will not discuss this situation in detail since the phase spaces we can find are analogous to 
the 1/ — 1 ones. 

The case of energy outflow is analogous to the small energy density analysis. In fact, the 
fixed point solution (62) can't be characterized by positive energy exchange parameter T. 
We can nevertheless have a critical point with energy outflowing from the brane into the 
bulk and negative Hubble parameter (as we can deduce from the expansion contraction, 
influx outflow symmetry). As the differential equation for p (59) shows, the energy density 
decreases and go to the trivial fixed point, meaning that the negative H critical point isn't 
an attr actor. 

The 4D energy density in the static compact extra dimension case is just given by a 
constant rescaling of the 6D density. Thus, the phase portraits are given by the plots in 
figure 1 for small energy density, and figure 2 for generic density (up to constant rescaling). 
However, for equal scale factors o? = 6, the 4D effective energy density q is dynamically 
determined by the energy density localized on the 5-brane p and the volume of the compact 
space V(2): Q — V{2)p. The 4D mirage density can be similarly defined as ;t = ^(2)X- We 
can single out the volume time dependence defining the dimensionful constant v such that 
V(2) = vb'^{t), where b{t) is as usual the compact space scale factor — in this particular case 
b{t) = a{t). The generic set of equations for the dimensionless variables H, g, i is 



where = A/v'^~^. We first note that in the case of zero energy exchange T = (j4 = 0) 
the 4D mirage density satisfies the 4D free radiation equation, as in the static internal space 
hypothesis. The 4D energy density q does not have a definite behavior in the case of energy 
outfiow. While in the static 2D compact space background it is clear that q, just as p, is 



^2 



Cp^d=6«^ _^ cy,d=6 (g + x) 



— K 



(80) 




(81) 
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suppressed in time since Wa=4 > 0, here we may have a negative coefficient for the linear 
term in g in (81). If < —3(1 + w)/2 — which is possible only for w < —1/3 if > — 1 
— we could have non trivial stable critical points, as in the energy influx context. This 
scenario would need further investigations. 

We can make some considerations regarding the more generic assumption on the relation 
between the Hubble parameters H and F, F — ^H, of subsection HID. For positive ^ 
the qualitative behavior is analogous to what we deduced in the case of static compact 
extra dimensions and equal scale factors. When ^ is negative, meaning that we are using 
a dynamical compactification approach [81], we could instead have some differences. In 
particular, it is worth noticing that (appearing in the conservation equation for p) can 
become negative, always implying a diverging behavior for the localized energy density at 
late time in the case of energy influx. Thus, there won't be stable critical point in the 
dynamical compactification scenario with energy flowing from the bulk onto the brane. 



C. Small density and free radiation equation of state: an explicit solution 

We can write an explicit solution to the set of equations (35)-(37) in the special limit 
of small localized energy density p <^ V and when Wd = d. This last condition is realized 
if -u; = 1/3 with static compact extra dimensions and if = (1 — 3w)/2 with equal scale 
factors. We must be carefuU though, because in this limit cv,d generally diverges. It is thus 
important to keep it finite by imposing a priori a specific value for w and w^^. With these 
assumptions, H"^ only depends on the sum (p + x) and the equation for this sum can be 
easily integrated independently of the explicit form of T. To be more specific, we get the 
following solution 

^ - (d-l)Mio ^^° + ^°^^"2(Z^6^ ^^^^ 
P + X = (po + Xo)^ (83) 

As in the four dimensional RS model with energy exchange analyzed in [43] , the evolution 
is determined by the initial value of the energy density (if we put, for example, xo — 0) 
weighted by the expansion in a (effective six or four dimensional) radiation dominated era. 

We deduce from (82) that \a{t)\ must have an upper limiting value for k > 0. For small 
positive a{t) the rate a{t) is a positive function and the universe expands until it reaches 
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the limiting value. If the compactification is over an hyperbolic space, the scale factor grows 
without bound. In particular, in a universe with extra dimensions evolving according to 
the same Hubble parameter as for the observed space-time, the expansion rate goes to a 
contant positive value as the scale factor grows. In a static extra dimension set-up instead, 
a{t) exponentially grows at infinity. When k = the explicit solutions for a{t) reduce to 
a{t) ~ i^/^ for static internal space, and a{t) ~ t^^^ for equal scale factors. These represent 
exactly a radiation dominated flat universe in four or six effective dimensions respectively. 
The evolution a{t) ~ t^^^ can also be traced in 4D to the Priedmann equation for scalar field 
subject to a null potential. 

If wc further assume the energy exchange to be linear in the localized energy density 
T — Ap, also imposing the initial condition Xo = and wa = d, the integration of the x ^^^^ 
p equations yields 

X = Po4(l-e-^% P = Po4e-^* (84) 

This solution shows that, for energy outflow A > 0, the initial amount of radiation energy 
density decays in favor of the mirage energy density. The late time evolution is thus governed 

by the mirage density. 

For equal scale factors it is interesting to write the explicit solution (remember Wd=6 = 6 in 
this case) in terms of the 4D densities q^tc {q — vb^{t)p, ;t = vb'^{t)x)- The ansatz implying 
static internal directions only results in a constant rescaling of the 6D quantities. The equal 
scale factor solution is given by 

^ - W^^^'^ + ^^^-To^ ^^^^ 

^ + ^ = (^0 + ^)4 (86) 

The evolution is still weighted by the characteristic effective 6D radiation dominated era 
But we can see, exploring the solutions for g and ^ in the case of energy exchange 
parameter determined by T = Ap (with positive A), that the 4D localized energy density 
evolves as 4D radiation, exponentially suppressed in time (as for the static internal space 
background). In fact, assuming for example to = 0, we get 




In this case A — A. 
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The case of energy influx is also analogous to the analysis in [43] . For both static compact 
extra dimensions and equal scale factors we rewrite the conservation equation for p for flat 
internal space as 

p + = (88) 

where we have used the equation determining H (82). This shows that for negative T, p 
should increase without bounds at late time. Since we are in the low density approximation, 
we can only rely on the generic p analysis of the previous section for large p. However, 
assuming T — Ap", we can deduce that for i/ > 3/2 the energy density can flow to zero (for 
certain values of the parameters). Indeed, u > 3/2 conrresponds to non stable crital point 
in the small density analysis. 

We remark that introducing the 4D density g for a universe with equal scale factors brings 

to 

Q + ^^Q^-T (89) 

Equation (89) tells us that the 4D localized energy density still grows unlimited at late 
time, if T is linear in g — more general considerations are analogous to the 6D density case. 
Again, we would need the full treatment for generic density. 
The acceleration parameter g = a/a in this context is equal to 

^^- (d-l)Mio (^0 + ^°)^ (90) 
For non zero k, the value of the acceleration can be either positive or negative. It has to be 
negative when k >0, but may be positive for compactification on hyperbolic spaces {k, < 0), 
giving as a result a loitering universe. 



V. CONSTRUCTION OF THE HOLOGRAPHIC DUAL 

The theory dual to the 7D RS model, via the AdS/CFT correspondence [58], will be 
derived in complete analogy to the 5D set-up considered in [10] . The RS model, with a time 
independent warped geometry, gives AdS^ metric as a solution to the equations of motion 
for the gravity action in the bulk. It will be useful to parametrize it according to Fefferman 
and Graham [90] 

GABdx^dx^ = jP'^dp'' + ep-^g^Ax^dx"" (91) 
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where the indices M,N run over the 7D bulk space-time, ji,!/, . . . span the 6D space- 
time on the 5-brane and p is a reparametrization of the z coordinate tranverse to the brane. 
The location of the brane, translated to this new set of coordinates, is p = which represents 
the boundary of the background (91). Generally, for all seven dimensional asymptotically 
AdS space-times the 6D metric g^j^^, can be expanded as [90] 

9 = 9(0) + P9(2) + P^^(4) + P^9(.6) + log p /i(6) + C'(p^) (92) 

The logarithmic piece appears only for space-times with an odd number of dimensions 
and is responsible for the cutoff dependent counterterm in the renormalizcd action. As a 
consequence, it is also responsible for the conformal anomaly of the holographic dual CFT 
[72, 76] — in fact, we don't have conformal anomaly in even dimensional backgrounds. The 
subindices in the coefficients of the metric expansion stand for the number of derivatives 
contained in each term. 

More precisely, RS background is a slice of AdSj, where the boundary gets replaced 
by the 5-brane and the IR part is reflected, eliminating the UV slice. To describe seven 
dimensional gravity we will take the usual Einstein-Hilbert action in the bulk, adding as 
usual a Gibbons-Hawking term [89] to take account of the boundary extrinsic curvature. 
Since the gravitational theory exhibits divergences in a space-time with boundaries, we 
also have to regularize the Einstein-Hilbert plus Gibbons-Hawking action, cutting off the 
boundary of the space-time. We are now going to illustrate the regularization procedure. 

A. Regularization on the gravitational side 

The renormalization for a gravitational theory in a background with boundaries has 
been explained in [72, 73, 76] for a generic number of dimensions. We will apply those 
computations to the case of a seven dimensional bulk space-time. 

In general, the bulk action for gravity gets modified with the Gibbons-Hawking boundary 
term [89] and with some counterterms also localized on the boundary 

Sgr — SeH + SqH — Scount (93) 

Using the Fefferman and Graham parametrization of the metric (91) and cutting off the 
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boundary ai p — e, we have 

Seh^M'J^ d'x^/^(^R[G] + ^y Sgh^2M'J^ d^xy/^K (94) 

where R [G] is the bulk Ricci scalar, K is the trace of the extrinsic curvature and jf^u is the 
induced metric on the boundary. Putting the brane at p = e corresponds to regularize the 
gravity action. The counterterm contributions necessary to make it finite in the limit e — > 
are given by 

S count — So + S1 + S2 + S3 (95) 

Si are terms of order i in the brane curvature R (the curvature of the induced metric 7^,^ 
on the boundary). In fact they can be written in terms of the induced metric 7^,^ and 
its Riemann tensor Rf^upa, using the perturbative expansion relating 7^,^ to g(o)^i, (see for 
instance [76]) 

So ^10 — J d6x^/=^ (96) 
S^^ --MH f d^xy^R (97) 

4 J p=€ 



s 



2 



^M'f ^_ d'xy/^ (r,,R^^ - ^R'^ (98) 



Ss = ^M5£5 / dSy=^ (IrR^^R'^'^ + + Ri^-RP-R^^^^ 



64 Jp^, " '\2 '^'^ 50 

+^i?'^'^V^V.i? - ^R^^'DR^,^ (99) 

The third order term S3 depends on the cutoff e and is thus responsible for the breaking 
of the scale invariance, i.e. it gives rise to the conformal anomaly for the dual 6D CFT in 
the context of the AdS/CFT correspondence. We also note that the zeroth order term is 
related to the brane tension term of the RS model (1) S^ens by Stens — ~2So, if we fine-tune 
^RS — 0. In fact, in the pure RS set-up, where the effective cosmological constant is null 
^RS — 0, the brane tension is V — 20M^/£, since the bulk cosmological constant is given 
by A7 = — 30M^/£^ as a function of the background length scale £ and the bulk Planck 
mass M, in our background metric parametrization (91). We will now use the AdS/CFT 
correspondence to compute the dual theory. 
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B. Gauge/gravity duality 



The AdS/CFT duality [58] is realized between gravity (string theory or M theory in the 
decoupling limit) in a background with one or more stacks of some kind of branes and the 
gauge theory that lives on the boundary of the near horizon geometry inferred by the branes 
[58, 59]. In our particular case, AdS7 x S"^ is the near horizon geometry of a system of 
parallel M5-branes in eleven dimensional M theory. The radius of the AdS space is given in 
terms of the eleven dimensional Planck lenght £pi and of the number A^ of M5-branes 

e = 2{T^Nfl^£pi (100) 

The radius of the four sphere is half the radius of AdS-j. The supcrgravity approximation for 
M theory is valid if A^ ^ 1 and ipi ^ N~^^^ — > 0, keeping the radius of the AdS large and 
finite in units oi £pi. The six dimensional theory that Maldacena [58] conjectured to be dual 
to M theory in the background described above is a (0,2) SCFT. This theory is realized as 
the open string theory in the world-volume of the M5-branes, in the low energy decoupling 
limit, and it does not contain dimensionless nor dimensionfuU parameters. The AdSr x 
supergravity backgroud is characterized by a 4-form flux quantized in terms of the number 
of M5-branes and is not conformally fiat, since the radii of the four-sphere and the AdS 
space are not coincident. 

The AdS/CFT correspondence relates the gravity (M theory) partition function for the 
bulk fields $j (which is a function of the value of the fields on the boundary of AdS^, (f)i) to 
the generating functional of correlation functions of the dual CFT operators with sources (f)i 

Zgr = J V^i e-^""- = Q-WcFT{<t>i) (101) 

Knowing that gravity on AdSr (the S'^ geometry can be factored out) corresponds to the 
specific CFT suggested by Maldacena [58], we can now obtain as a consequence the theory 
dual to the 7D RS model, in analogy to [10]. In fact, the action of the gravitational theory 
that we want to analyze via holography is 

SrS — SeH + Sgh + Stens + (102) 

We just add the Gibbons-Hawking term to (1). We expect that the hidden sector of the 
holographic theory reflects the bulk non trivial contents encoded in the bulk components of 
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the 7D stress-energy tensor (Tj", Ty ) when we go to the non conformal interacting general- 
ization. 

We now have to keep in mind that the duahty for gravity on AdS^ can be stated as 

Jp>e 

Secondly, we have to remember that Stens — — 2S'o. Furthermore, we note that the integration 
in (103) is over one half of the space-time appearing in the RS model, because of the 
Z2 reflection along the z direction. Since the integrals over the two specular regions are 
independent and equal we can write 



J alio 



'allp 
Jp>e 

where Xi ^^e the matter fields on the brane. Finally, putting all toghether, using equation 
(103), we obtain 



Zrs [(pi, Xi] = [ 'D^iVxi e-^^cFT-2S,-2S,-2S^-s^ (^q^^ 
The RS dual theory is 

SrS = ^CFT + Sr + Sr2 + Sr3 + Sm (106) 

having defined 

ScFT — "^WcFT, Sr — 2Si, Sr2 — 2S2, Sr3 — 2Ss (107) 

The 6D Planck mass is thus given by Mpi — 

We are now ready to calculate the equations of motion for the holographic 6D RS cos- 
mology. 

VI. HOLOGRAPHIC EVOLUTION EQUATIONS 

As we know, the RS classical solution in a 7D bulk with a warped geometry is AdS7. Since 
our purpose is to study the cosmology associated to the 7D RS set-up, we have generahzed 
the ansatz for the metric to be time dependent in section II and III. We have successively 
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reviewed the notion of holographic dual theory in the previous section. What we want to do 
now is to describe the cosmology of the seven dimensional RS model from the six dimensional 
holographic point of view, using the correspondence relation obtained in the previous section 
and generalizing the ansatz for the 6D induced metric on the 5-brane to a time dependent 
geometry, as we had done for the 7D bulk analysis. 

We consider a 6D space-time, compactified on a 2D internal space, with a FRW metric 
for the four large dimensions. The induced metric tensor can be expressed as 

7^^dx''dx^ = -dt^ + r^^x^dr^ + a\t) r^d9^ + a\t) sin^ Odcj)^ + 

J. thV 



^T^^.dp' + b\t) pM^' (108) 

1 — Kp^ 



where k and k, a{t) and b{t), H{t) and F{t) are, respectively, the curvatures, the scale 
factors, the Hubble parameters for the 3D and 2D spaces. 
The action we are considering is 

^RS = ScFT + Sm + Sx + Sii + 5*^2 + Sr3 (109) 

Sr, Sr2, Sr3 being respectively twice the first, second and third order terms in the curvature 
contributing to the counterterm action, as defined in (96)-(99). Sx is an effective cosmo- 
logical term on the branc — that represents a generalization to the case of a non exact RS 
fine-tuning with respect to the action (106). and Scft are the matter and (twice the) 
CFT actions of the 6D description. 

We want to solve the Priedmann equations, imposing the conservation and anomaly 
equations, defining 

T,. = -L^ W,^ = 1 ^ (110) 

and V^j,y = W^y + Y^^ + Z^y. The equations of motion take the form 



= (112) 
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Here A(6) is the general anomaly for a 6D conformal theory [75, 78] that comes uniquely 
from the 5*3 contribution to the renormahzed action [72] -[76] [94], while Y is the trace Y^f 
of the variation of (twice) the second order counterterm action Sji2. The trace of Z/^^, is 
null [95]. The trace of Yj^i, is quadratic in the curvature of the metric (108) 

Y = ^M'i (^R^'^R,, - ^R'^ (113) 

The explicit form for the anomaly is a complicated expression of dimensions 6, cubic in the 
curvature, and is discussed in appendix A. The effective cosmological constant on the brane 

is A. The stress-energy tensors are parametrized as 

700 = P(^), Tij ^ p{t)^ij, Tab ^ T^{t)lah ^^^^^ 
Vqq = (T(t), Vij = Upit) 7ij, Vah = (T^{t) Jab 

where the indices ij... parametrize the space part of the 4D FRW space-time and run from 
1 to 3, while ab . . . belong to the 2D internal space and take values in (4, 5). [96] 

Equations (112) take the following form when we choose the metric (108) and the stress- 
energy tensors written in (114). The Priedmann equations become 

Mj^i (sH'' + 6H F + F'' + 3^ + =p + a + X 

(2H + 3H'' + AH F + 2F + 3F^ + ^ + = -p-ap + X (115) 

M^; (sH + 6H'^ + 3HF + F + F^ + ^^ = -tt - cr^ + A 

then the conservation equations 

& + 3{a + ap)H + 2{a + a^)F = 

P + 3{p + p)H + 2{p + tt)F = (116) 

and finally the anomaly equation 

a-3ap-2a^ = A(Q)+Y (117) 

As we said, the anomaly comes from the cubic counterterm, so that it is cubic in the 
curvature. We make some more precise statement about its form in appendix A, where we 
also explicitely give Y. 
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A. Simplifications and ansatz 



The set of equations (115)-(117) does not contains six independent equations. Plugging 
the conservation in the first Priedmann equation differentiated w.r.t. time, we get a hnear 
combination of the other two Priedmann equations. So we will discard the last of (115) 
from now on. We further note that the system contains only one algebraic equation: the 
Priedmann equation. We will start by solving the anomaly equation in terms of one of the 
pressures coming from the hidden theory. 

Plugging the expression for Up obtained evaluating (117) into the first of the conservation 
equations (116), we get a differential equation for a depending on cr^ 

(7 + 2 (2// + F) a + 2 (ii- - F) (7^ - >1(6) + r (118) 

To obtain a solvable decoupled equation for u we can consider the limit in which the internal 
space has the same CPT pressure Ot, = ap of the three large dimensions [97]. Else, we can 
also consider the limit of zero pressure — for the CPT — in the internal space. Putting 
these two limits toghether, we can try to solve the Priedmann equations imposing a more 
general ansatz 

(Ttt = flap (119) 

So that 

(7 - -(7p = ^(6) + Y (120) 

where cu = 1/(3 + 2^2) (a; is equal to 1/5, 1/3 in the two limits considered above) and the 
differential equation for a becomes 

(7 + 3 [(1 + uj)H + (1 - uj)F] a-[3ujH+{l- 3uj)F] (^(g) + F) = (121) 

We could now evaluate a solving the following integral 

a3(l+u))53(l-(^) J 

• [caF(6)+Cb/(6) + F] (122) 

where x is a solution for the homogeneous equation 

X + 3[(l+c.)// + (l-c.)F]x = X - ,3(i4b3(i-.) (123) 



3lu- + 1 - 3a; - 
a b 
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We observe that generally (122) is not explicitely integrable. In (122) we have written the 
anomaly in terms of its contributions that are the Euler density in six dimensions £^(6) and 
the local covariants included in /(g) (see appendix A for further details). 

The set of independent equations we finally have to solve, once we use (120) to eliminate 
Up by means of 



-cu(^(6)+F) (124) 



is then 



M^;(^3i/' + 6//F + F' + 3^ + p^ =p + (7 + A (125) 

M^Pi (2H + 2F + 3H^ + 4HF + 3F^ + ^ + = -wp - cua + uj {A^e^ + Y) + X 
(j + 3 [(1 + uj)H + (1 - u)F] a = [3ujH + (1 - 3u)F] (^(g) + Y) 
p+[3{l + w)H + 2{l + w^)F]p = (126) 

(the anomaly ^(g) will be written explicitely — in terms of H and F — in the particular 
cases that we will take under examination in the following). We also use the three following 
ansatz relating the pressures and the energy densities 

p = wp 

t: ^Wp^w^p (127) 
a^^VtUp (3 + 2Q = l/a;) 

Now we are left with a system of four equations (125)-(126) in four variables {H,F, p,a). 
The other variables (cTp, (TttjP, tt) are determined by the ansatz (127) and by the equation 
(124). In the next section, this system of differential equations will be studied restricting to 
some special limits, such as flat or static internal space, or equal scale factors. We will find 
the critical point solutions and analyze the associated stability matrix. 



VII. HOLOGRAPHIC CRITICAL POINT ANALYSIS 

The fixed points of the cosmological evolution of the universe we are considering may 
represent its infiationary eras — for instance the early time or the late time acceleration 
— , since the Hubble parameters, just as the energy densities, are constant. If the constant 

43 



value for the Hubble parameter is positive we have inflation. In this section, we are going 
to look for the existence of such inflationary points for our specific holographic model and 
to find what kind of dependence they have on the parameters of the theory. 

We will describe the fixed point solutions in the special limits of flat extra dimensions, 
curved static extra dimensions and equal scale factors for the internal and extended spaces. 
We will then study the stability matrix associated with the critical points. Since the flxed 
points represents inflationary eras in the universe evolution, they could offer an explanation 
to the early inflation or to the late time acceleration. In the flrst case they will have to be 
unstable or saddle points to allow the trajectory of the cosmological evolution to flow away 
from inflation and exit from this phase. In the second case, on the other hand, the fixed 
points must be stable and act as attractors for the nearby trajectories. 

In what follows we will always suppose that the effective cosmological constant on the 
brane A will be zero, unless we specify it differently. 



A. Flat compact extra dimensions 

In the limit of zero spatial curvature for the extra dimensions and for the extended space, 
the Friedmann plus conservation equations (125) take the form 

Mj^i {3H^ + 6HF + F^) ^ p + a + X (128) 

Mj>i (2H + 2F + 3H^ + AHF + SF^) = -wp - uja + uj {A(&) + Y) + \ (129) 

(7 + 3 [(1 + u)H + (1 - u)F] a = [?,uH + (1 - 3u;)F] (^(g) + Y) (130) 

p+ [3(1 + «;)// + 2(1 + «;^)F]p = (131) 

We find the fixed points of this system of differential equations and study stabihty with 
some further restrictions (see appendix B in the flat extra dimension subsections for the 
explicit calculations). As we point out in appendix A, the anomaly A{(o) generally depends 
on the Hubble parameters of the model, on their time derivatives up to the third order and 
on the spatial curvatures. This remains true also for the flat extra dimension limit that we 
are examining, ignoring the curvatures. 

Fixed point solutions With the assumption of flat internal space and zero curvature for 
the 3D space as well, we can flnd different flxed points depending on the value of the extra 
dimensions Hubble parameter. They can be summarized as follows. 
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1. As we illustrate in appendix B, there are two non trivial time independent solution with 
F non vanishing at the fixed point, F^, ^ 0, one for a; 7^ 1/5 and one for a; = 1/5. The 
values of the 3D Hubble parameter (the measurable Hubble parameter) are given in 
terms of the constants cj, c^, c^, cy and the mass scale Mpi. We note that the anomaly 
parameters c^, Cb are given by the CFT, while u relates the hidden sector pressure of 
the internal space to the hidden pressure of the 3D space (127). In particular, the two 
Hubble parameters are related by the equahty if* = (Cg + 1)F* (where is a rather 
complicated function of u!,ca,cb,cyMpi) when cu ^ 1/5, while for a; = 1/5 we have 
Hi, — Fi,. For (jj — 1/5, Fi, ^ and we can choose it to be positive or negative, implying 
that the extra dimensions scale exponentially at those points, with respectively either 
positive or negative velocity. Consequently H would also describe either a contracting 
or expanding universe. 

2. We also found a fixed point solution for which the extra dimensions are static, i.e. F* = 
(a; 7^ 0, —1 and cb 7^ 0), meaning that, while our visible universe is exponentially 
growing (or, in principle, decreasing) the internal space isn't expanding nor collapsing. 
The corresponding solution is given by 



20 , ^9 

Hi = M% 



48cy ± W6 ( 3844-cb ^ 



UJ+1 



(132) 



2lu 20 00 

-O-TT* = -CTn* = rr^Pl 

3a; — 1 cbuj + 1 





48cy ± ( 3844-cs ^ 



u;+l 



The roots are real if 384cy — cboj /{uJ + 1) > and cannot be both positive. We never 
have a couple of fixed points in the phase space diagram. 

3. A third fixed point is characterized by zero extra dimension Hubble parameter and 
Lo — —1. In this case the critical point exists only if the conformal field theory is 
characterized by a positive coefficient for the type B anomaly. The solution is 

^,^ 640cy a. = -a,. = i<7.. = ^M|„ p. = (133) 

Cb 2 3cb 

If Cb is zero (i.e. the anomaly vanishes at the fixed point) we are left only with the 
trivial fixed point. 
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4. For vanishing a;, A should be non zero to get the inflationary fixed point — 

5. There also exists a trivial fixed point with = F^^ — 0, where also the anomaly and 
the trace become zero, and p^, — a^, — if cu ^ w or — — cr* if cu — w. 

In any case — i.e. for every uj and A — the solution docs not depend on w, (in fact, the 
system of equations (B1)-(B5) doesn't contain w, w^). So, if a solution exists for some cu and 
ca,cb,cy, that solution always exists whatever values the two parameters relating matter 
pressures to energy density take. This marks a difference with the bulk analysis of section 
IV, since here we don't have any bulk dynamics perturbing the conservation equations, being 
the hidden sector theory conformal and non interacting. 

All the critical points have zero localized energy density p^ (except for the trivial point 
with u) ~ w). Also, when the Hubble parameter is non vanishing, we don't get a positive 
value for p^. 

Stability analysis In appendix B we analyze the stability of the = critical points 
linearizing the system of differential equations around the fixed point. We conclude that 
the studied fixed point (132)-(133), characterized by vanishing and a; = —1 or a; 7^ —1 
can both be saddles or attractors, depending on the value of the anomaly parameter Cb, of 
Cy and of a; (relating the two CFT pressures). The trivial fixed point cannot be analyzed 
at linear order, since its stability matrix is null. It is obvious from (131) that for positive 
Hubble parameters the energy density goes to zero at late time. 

We can thus observe that, starting the cosmological evolution with a hidden energy 
density a different from zero and with a suitable value of the anomaly coefficient cb, choosing 
UJ to be such that the = stability matrix has negative eigenvalues, the corresponding 
F^ = fixed point could be a global attractor for the fiat extra dimension universe. This 
critical point could eventually represent the present accelerated era. However, it's a zero 
energy density critical point. 
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B. Static compact extra dimensions 



The Einstein equations of motion (125) get simplified when we take the b — const ansatz 
for the internal space scale factor 

(134) 

(135) 
(136) 

(137) 



/)2 



(2H + ^) = -'^p - + ^ (As) + y) +^ 

a + 3(1 + u)Ha = (^g) + Y) 

p + 3{l + w)Hp = 



Prom these equations we can deduce the corresponding fixed points and their criticality, 
following the calculations in appendix B. 

Fixed point solutions We have found the inflationary fixed points for a universe with non 
evolving internal dimensions, i.e. with constant scale factor b{t) = bg. Besides the trivial 
fixed point — there are other solutions. 

1. The existence of non trivial fixed points is determined by the values of the parameter 
of the specific conformal theory ca,cb,cy and lu, but also by the mass scales Mpi and 
K/bg. The form of the fixed point solutions are derived in appendix B. An easy critical 
point solution can be derived in the case of zero type B contribution to the anomaly 
and cu — —1 



tt2 

= I2 
^0 



—(7, 



9± W90 



5ca 1 K, 
192cy 6| 



(138) 



p* 



I 3 



9± 4 90 + 



5ca 1 K 
192cy 6g 



K 
1? 







For K = it reduces to the trivial fixed point. We have real roots for if 192cyMp; > 
— c^/t/456g and they will be both positive when k > and —192cyMp^ < CAn/4:5bl < 
— 192cy M PI / 252, so that ca must be negative. 

2. We also have a trivial critical point with = 0. In particular both p^ and cr* can 
be different from zero, they are functions of n/b^ as we can see from the equation 



47 



p^-\- a-i, — MpiK,/bl, and are also related by 

(l^„),.^(l+.K.„(_i^| + !|:M|,)| (139) 

3. Other fixed points can be found, for example for cu — 0. The internal space curvature 
K would have to be negative in those cases, since — —K/Zb^ and we would have to 
compactify on an hyperbolic space. 

Stability analysis For the static extra dimension fixed points we found that, depending 
on the values of c^, c^, cYMpi and a;, we can get an attractor or a saddle. 

We can thus choose the hidden sector parameters such that we can get a stable fixed 
point. There exists however another critical point, i.e. the trivial one characterized by 
Hi, = Q but generally non zero and a*, which always is a saddle. So, trajectories may 
either be attracted by the non trivial critical point or flow away from the saddle point. 



C. Equal scale factors 

Another limit that simplifies some of the calculations is the equal scale factor assumption. 
In this case the Hubble parameters of the internal space and the 3D space are equal, F = H, 
and we remain with the following set of equations for the variables H, p, a 

M^; (lOH^ + = p + a + X (140) 

M^i (ah + lOH^ + ^) = -coa + LU {A^^^ + Y)+X (141) 

a + 6Ha = H (^(g) + Y) (142) 

p+ {5 + 3w + 2w^)Hp = (143) 

Fixed point solutions We find an infiationary fixed point where the internal space is also 
staying in an infiationary era, since the two Hubble parameters are equal (see appendix B 
for calculations). The existence of such fixed point depends on the values of ca,cb,cy, h 
and uj. 

Is is necessary to impose uj — 1/5, implying that the two pressures characterising the 
CFT stress-energy tensor must be equal (Tt^ — Up, in order to obtain the following fixed 
point solutions. 
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1. We write the explicit solution for a flat internal space k, — 0, where the critical value 
for the Hubble parameter is given in terms of the anomaly coefficients and of the 
only dimensionful parameter which is the 6D Planck mass. This solution restricts the 
possible values of ca_ and cb, coming from the conformal field theory anomaly 

24 



ca + 2cb 



24cy ± 
240 



5764 + (CA + 2cb) 



-Ml,, 



24cy ± J5764 + {CA + 2cb) 



(144) 



ca + 2cb 

For (24cy)^ > —{ca + 2cb) the two roots are real. We must satisfy the condition 
{ca + 2cb) < on the anomaly parameters in order to have two positive if* critical 
points (both denoted by zero energy density Pi,). 

2. For K, there exists no fixed points, since (.4(6) + and [lOH^ + ^) cannot 
be both constant unless we enforce the staticity condition on the Hubble parameter 
= and K = 0. 



Stability analysis In this last case of equal scale factors we have again carried the stability 
analysis in appendix B. As a result, we found that the equal scale factor critical point with 
flat internal space is an attractor. It could thus represent the eternal acceleration of the 
universe. 



D. Comments 

We want to summarize the interesting features of the critical analysis of sections VII A- 
VII C, in the perspective of the comparison with the 7D bulk gravitational dual description 
that we illustrated basically in sections HI and IV. 

- All the critical points we can find in the brane description are characterized by ex- 
actly zero value for the localized matter energy density p (except for some H^, — 
trivial points). To have a non vanishing energy density it is necessary to introduce an 
interaction term between the matter fields and the hidden sector fields. The reason 
for this is that it modifies the conservation equation for p, allowing for a non zero 
time independent solution. This intuitively corresponds to turning on the brane-bulk 
energy exchange on the bulk gravity side. 
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- In most of the simple critical point solutions, the hidden sector pressure api, is related to 
the energy density a^, by — — cr*. This indicates a vacuum behavior for the equation 
of state of the hidden sector of the holographic dual theory at the inflationary fixed 
points. 

- There also appears to exist more than one critical point solution in some of the ex- 
plicitely examined limits. Since the stability matrix analysis reveals that we can have 
either stable or saddle points (depending on the CFT and counterterm parameters 
ca,cb,Cy and on the Plack mass), we can expect two kinds of behavior (if the number 
of critical point solution is two). Whenever one of the fixed points is attractive and 
the other one is a saddle, we can generally depict a phase portrait such that some of 
the trajectories are attracted by the stable point, while others can be repulsed by the 
saddle and go toward the large density region. If, on the other hand, we get two saddle 
points, trajectories bend near to the saddles and flow away. We note that the trivial 
critical point has undefined stability at linear order in the perturbations, so that it 
may either attract or repel trajectories in its neighborhood. However, ii H > and 
WjWt^ > —1, late time evolution always is described by p ^ 0. 

- Comparing these results with the bulk cosmology we note that we get both trivial and 
non trivial accelerating critical points in the CFT description. The non trivial fixed 
points are associated to the so called Starobinski branch of the solution to the Ein- 
stein, conservation and anomaly equations. Acceleration is produced due to conformal 
anomaly in higher derivative terms. These accelerating points can be interpreted to be 
in correspondence with the accelerating fixed points found in the gravity description 
with infiow. The trivial fixed points come from the smooth branch, instead, and they 
are smoothly connected to the gravity description. An analogous map can be found 
in the 5D-4D model [10]. 

- In the bulk description we only found positive acceleration critical points since q = H^. 
With the holographic approach, we could in principle also get non constant H critical 
point solutions. In fact we could solve the system of Einstein equations on the brane 
asking time independence ior q — H + H'^,p and a. We would get a new system of 
first order (non linear) differential equations in H, which in principle could have non 
trivial solutions. 
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VIII. BRANE/BULK CORRESPONDENCE 



We derive some solutions that can illustrate interesting aspects of the cosmological model 
we considered and in particular of the duality that relates the two descriptions. We will be 
able to find explicit expressions by making special simplifying assumptions on the parameters 
of the holographic theory and on the space-time background. These examples allow us to 
make a comparison between the results we will find in the holographic set-up and the 
expressions we derived in the bulk gravity theory. 

Since we are interested in comparing the two dual approaches, in the sense of AdS/CFT 
correspondence, we have to derive some expressions for H in terms of the localized matter 
energy density p and of a mirage density x- In section III, we performed the cosmological 
analysis in the 7D bulk description, expressing the 3D Hubble parameter H in terms of the 
localised energy density p, reducing the first order ODE in H to an algebraic equation for 
plus a first order ODE for the mirage density %. In the holographic dual theory, the 
mirage density is identified with the solution to the homogeneous equation associated to 
the conservation equation for the hidden sector density a. It will thus have the property of 
obeying to the free radiation conservation equation in d effective dimensions (where d is the 
effective number of dimensions equal to 4 in a static compact space background and to 6 
when the a{t) and b(t) scale factors are equal). 

To obtain the explicit result for in the brane dual description, however, it is necessary 
to integrate the differential equation (122) for the energy density a. It is established that 
the anomaly and the trace of the quadratic contribution to the variation of the dual theory 
action are highly non trivial functions of the Hubble parameters and the spatial curvatures 
and contain derivatives of H, F up to order three. So, an analytical integration of the 
a conservation equation is in general apparently unachievable. However, it is possible to 
neglect the ^(e) and Y contributions if we are in the slowly scaling approximation, which 
corresponds to a small curvature approximation. This is what we are going to discuss in the 
following. 
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A. Slowly scaling approximation 

We give a rough idea on how the correspondence between the brane and the buUk dual 
theories works. In fact, we will neglect all the higher order terms in the holographic descrip- 
tion, which is equivalent to ask that the Hubble parameter is negligeable with respect to 
the Planck mass <C Mp^ (i.e. a -C Mpio). In this approximation, all the higher order 
curvature terms — including the anomaly and the trace Y"^ — can be neglected in favor of 
contributions proportional to the Einstein tensor. The integration of the a equation would 
give terms of the order of MpiH'^ (from Y) and (from A(6)). Once we plug the result 
for a in the Einstein equation (125) these terms are suppressed, since the l.h.s. of (125) is 
of order M^;//^ 

We anticipate that, as a consequence of the small Hubble parameter approximation on 
the brane, we get a linear dependence of on the mirage plus visible matter energy 
densities (the hidden sector density a is identified with the mirage density x) . Higher order 
contributions due to the anomaly and to Y would give rise to higher power dependences in 
a small density expansion for H^. Since in the holographic description we truncate to linear 
order in density, we also keep only linear terms in p for the bulk gravity results. The bulk 
equations for H and x (35)-(37), can be formulated independently of w if we ignore higher 
(quadratic) order terms in p and assign a specific value to -u;^. Neglecting the second order 
term in p we will only have one condition to determine M and V in terms of the brane 
parameter Mpi, so that only the ratio M^°/V will be identified. 

Since in this approximation the quadratic and higher order dependence of H'^ on p are 
absent, we don't capture the eventual Starobinsky [79] behavior of the solutions to the 
Einstein equations [80]. The higher derivative terms are necessary in that case to calculate 
the exit from inflation to a matter dominated universe and the subsequent thermalization to 
radiation dominated era. In Starobinsky model the higher derivative terms are represented 
by the type D conformal anomaly contribution to the trace of the stress-energy tensor. 
Nonetheless, in the 5D RS holographic dual analysis of [10], where these terms are cancelled, 
stringy corrections like Gauss-Bonnet terms can play the same role. We note that in our 
set-up, the 6D conformal anomaly contains suitable higher derivative terms, not only in the 
total derivative contributions but also in type B anomaly. 

We now explain the results to which the slowly scaling approximation leads in some 
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particular limits. 



1. Equal scale factors 

We start looking at the Einstein, conservation and anomaly equations in the equal scale 
factor limit. The system of equations (125)-(126) for the theory on the brane takes the form 

Mt>i (lOH^ + 3^ + ^] = p + a + X (145) 

Ml.1 [AH + lOH^ + 4 + 4 ) = - + ^ (Ae) + ^) + A (146) 

& + 6Ha = H (^(6) + Y) (147) 
p+(3(l + w) + 2(l + w^))Hp^0 (148) 

We note that the set of equations is independent of the hidden sector parameter cu and 
it is furthermore interesting that the homogeneous equation associated to (147) is indeed 
precisely the 6D free radiation equation, independently of the value for cu. In particular, the 
two conservation equations can be written in the integral form 

a^X + a-' I dta'H{A(e)+Y), x = Xo (^)' (149) 



3(l+«;)+2(l+«,^) 

a J 



Plugging the result for a in the first equation of the system (145) and neglecting the curvature 
higher order terms that come from the integration of (^(e) + ^) > "^^ obtain the following 
expression for H^, toghether with the p and x equations in their differential form 

X + 6Hx^0 (151) 

p+{3{l + w) + 2{l + w^))Hp = (152) 

It is now easy to compare (150)-(152) with the corresponding system of equations in the bulk 
description of the equal scale factor universe, with zero brane-bulk energy exchange Tj and 
bulk "self interaction" Tj. The expression for can also be written in a w-independent 
way fixing ty^ (in this particular case we could also include the quadratic term in p in the 
w-independent formulation, i.e. when — w). Neglecting the second order term in the 
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energy densities we obtain 

X + 6Hx^0 (154) 
p+ (3(l + w) + 2(l + w^))i7p = (155) 

The two systems of equations (150)-(152) and (153)-(155) perfectly agree at this order 
in the approximation. The matching between the scales on the two sides of the duality is 
then 

_ /IS /\/f4 — PI n 

-y- - ^cv(,q)Mpi — > - [iob) 

As we announced, only the ratio M^^/V can be determined, sice we only have one condition 
to match the two descriptions. When higher order corrections are included in the brane 
description wc would generally find a matching for both M and V, which in principle would 
depend on the particular CFT parameters {ca, cb, cy). We can guess that the ratio M^^/V 
would not depend on them (indeed, for the pure RS set-up we get MpjV oc M^°). When 
— w (i.e. when the pressures of the matter perfect fluid relative to the 2D internal space 
and the 3D space are equal tt — p) the coefficient cy becomes cv(eq) — 1/80. It is interesting 
to note that for w = w-^ the matching exactly reduces to the RS condition for zero effective 
cosmological constant on the brane Xrs = 0. Since we are in the limit F = H, it seems 
natural to have tt — p too. 

2. Static compact extra dimensions 

We consider the static extra dimension hmit F = 0. The Einstein equations plus conser- 
vation and anomaly equations in this limit read 

k K 

^0 



Ml.i[3H' + 3— + -^]=p + a + X (157) 

(2H + 3H^ + \ + ^) = -wp-u;a + u{A(^e)+Y)+X (158) 

a + 3(1 + uj)Ha = ScuH {A(6) + Y) (159) 

p + 3{l + w)Hp^0 (160) 



54 



The homogeneous equation associated to (159) is the 4D free radiation equation only if 
cu — 1/3, implying that the hidden sector pressure of the internal space (7,r must be zero. 
With this assumption, the two conservation equations become 

a^X + a-'J dta''H{A^e) + Y), x = Xo (^)' (161) 

/ao\3{i+"') 

^ = u 

The results for x ^-nd p agree with the bulk formulation for zero energy exchange. Plugging 
(161) into (157) and neglecting the curvature higher order term, as we are in the slowly 
scaling approximation, we find 

X + ^Hx^O (163) 

p + 3{l + w)Hp = (164) 

Although conservation equations agree, the Priedmann like equation doesn't give the ex- 
pected 1/6 coefficient in front of the n/bl term in (162) that instead follows from the equa- 
tions on the bulk gravity side 

X + AHx^O (166) 

p + 3{w + l)Hp = (167) 

The bulk equations are derived in the density linear approximation and for vanishing energy 
exchange T° and T/ . The coefficient cvist) can be written in a w-independent way if we fix 

As a consequence, the ratio of the two bulk parameters can be identified with 

The n/bl terms differ in the two dual descriptions (in the static extra dimension background). 
The matching (168) gives a result that depends on the valus of w^w-,^ in a different way if 
compared to the equal scale factor limit (156). It is thus interesting to further examine 
how the matching varies according to the value of the internal space Hubble parameter. We 
are indeed going to consider the proportionality ansatz F — to better understand this 
behavior. We note that in the hmit w^^ — {w + 5)/6 we recover in (168) the RS fine-tuning 
determining zero effective cosmological constant on the brane Xrs — 0, since cv(st) = 1/40. 
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3. Proportional Hubble parameters 



Following the computations in the last two sections and generalizing them, we derive 
the set of equations for H'^,x and p for proportional and small Hubble parameters. Since, 
as before, cr = x if we neglect higher order terms in the small curvature approximation, 
equations (125)-(126) lead to the following equations 

1 Kb 1 

^ ie^ + ee. = ie^ + 66 + 3)M|, + ^^^^^ 

X + d^,Hx^O (170) 
p + w^^Hp = (171) 

where = 3(l+a;) + 3^(1 —cu), w^^^ = 3{l + w) + 26(l + u'7r) and ^ is the proportionality 
factor F = (or b = a^''). The bulk equations where derived in (38)-(40). For the moment 
we will keep two different proportionality factors: in the bulk F = ^bH as in section III D 
(putting zero T/ and Ty ) 

+ (el + 3^B + 6)^= (2eB + 3)Mio(^ + ^) (1^2) 
X + 4b^X = (173) 

p + w^^Hp^Q (174) 

We have used the following definitions: ci^^ = 6(,^| + 2,^^ + 2)/{2^b + 3), = 3(1 +w) + 
2'CB(1 + W7r) and we recall that cy_^ is cy.j = cv/{w^B — d^B), where cy = (31^1 — Gwtt — 5)/400. 

In order for the two descriptions to be equivalent w.r.t. the x and p differential equations 
(which don't get any correction from higher order contributions), we have to put — Cb — 
as it was expected. Besides, the parameter relating cttt to dp must be a; = 1/(2^ + 3). 
However, assuming an equal proportionality relation on the two sides of the duality, the 
coefficient of the k term in (169) and (172) differ if the two curvatures in the bulk and brane 
descriptions are equal, unless ^ = 1. So, the only set-up that predicts the same effective 
spatial curvature for the internal space in the brane and bulk descriptions is the equal scale 
factor background (neglecting higher order corrections). However, we can determine an 
effective spatial curvature for the internal space in the brane description, given by Kb — 
(e' + 6e + 3)«B/(^' + 3e + 6). 

The matching for the scales of the two dual theories is given by 
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-2cv,,Mt.i (175) 



It is always possible to choose a such that the matching relation (175) gives the RS 
fine-tuning condition Mpi — 20M^°/V. Otherwise, missing the fine-tuning would amount 
to introducing a non vanishing effective cosmological constant on the RS brane. 

IX. NON-CONFORMAL AND INTERACTING GENERALIZATION 

To examine the general cosmological evolution that reflects the presence of energy ex- 
change between the brane and the bulk in the seven dimensional bulk picture and the non 
zero value of the bulk component of the stress tensor Tj, we will drop the assumption of 
having a conformal non interacting field theory living on the brane. Intuitively, a non van- 
ishing T° in the bulk description corresponds to interactions between the gauge theory and 
the visible matter. The diagonal Tj component appears in the brane description as dual to 
a new trace term spoiling the conformal invariance. The generalization of the 6D RS dual 
action will be modified adding the new interaction term and substituing the hidden sector 
CFT with a strongly coupled gauge theory (SCOT). The following analysis will be done in 
analogy with the 4D holographic dual generalization of the 5D RS cosmology exposed in 
[10]. 

Using the notations of section V we write the generalized action as 

Sgen — SsCGT + Sji + Sff2 + Sffs + Sm + Sint (176) 

where the new entry is the interaction term S'j„t and Soft has been changed into Sscgt- 
The strongly coupled fields can be integrated out, transforming the sum of the strongly 
coupled theory action plus the interaction term into an efi^ectivc functional of the visible 
fields (and of the metric) Wscgt- As a result, the action (176) becomes 

Sgen = Wscgt + Sr + Sr2 + Sr3 + Sm (177) 

As in the conformal non interacting case, we are now ready to calculate the general 6D 
equations of motion for the holographic generalized RS cosmology. 
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A. Generalized evolution equations 



The stress-energy tensors are defined in an analogous way 

1 SSm 1 SWsCGT 



T — — W 



7 57^*^ <^7' 
1 5Sji2 1 SSrs 

'1X1/ ^' 1X1/ ^ ■'■ 1X1/ ~ ^ixi/t ^ ^ — 

and the Einstein equation, the (non) conservation conditions, the anomaly equation read 

^iv ' ■'■ ixv ' * ixv 





= T ■ 

-'- /xv 




= T 




= -T 




= ^(6) 



X + Y (179) 

The total stress-energy tensor is still conserved. Taking account of the interactions between 
the hidden theory and the matter generally amounts to have non separately conserved V^,^ 
and Tfj^^. This is reflected by the introduction of a non homogenous term in the conservation 
equations. The anomaly equation contains the general expression for the conformal anomaly 
in six dimensions A[6) and the trace term Y. Furthermore, it gets modified including an 
extra term X that accounts for classical and quantum breaking of the conformal symmetry 
in a FRW plus compact space background. The stress-energy tensors are parametrized as 
before 

Too ^ pit), Tij ^ p{t)jij, Tab ^ T^{t)lah _ 

Vbo = o-(^), Vij = ap{t)jij, Vab = cr^{t)jab 

The consequent changes in the equations written in terms of the Hubble parameters, of 
the energy densities and pressures are the following. The Priedmann equations remain the 
same (as for the non interacting conformal case) 

(sH^ + 6HF + F^ + 3^ + ^^ = p + a + X 

M^i (2H + 3H^ + 4HF + 2F + 3F^ + ^ + ^^ ^ -p-ap + X (181) 
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the conservation equations now involve the quantity T 

d + 3((T + ap)H + 2{a + a^)F = T 

p + 3{p + p)H + 2{p + 7r)F ^ -T (182) 

and the anomaly equation includes the conformal breaking term, as a consequence of the 
masses and /3-functions of the strongly coupled gauge theory 

(7 - 3(7p - 2(7^ = ^(6) + Y + X (183) 

X has to be written in terms of the /5-functions and operators of the SCGT and matter 
theory. Taking the same ansatz for the pressures as for the non interacting conformal theory 

= nap = 3 + 2Q) 

p = wp, TT = Wt^P (184) 

and using the anomaly equation to eliminate ap — cu (a — A^e) —Y — X^ from the set of 
remainig equations, we get 

Ml,i(3H^ + 6HF + F^ + 3^ + ^^ =p + a + X (185) 

Mj^i (2H + 3H^ + 4HF + 2F + 3F^ + ^ + ^^ ^ -wp -cua + X + uj (A^e) + Y + X) 

(186) 

a + [3(1 + u;)H + 3(1 - u;)F] a = [3u;H + (1 - 3u;)F] {A^e) +Y + X)+T (187) 
p + [3(1 + w)H + 2{1 - w^)F] p = -T (188) 

The cosmological evolution described by these differential equations which include non 
conformality (represented by the X term) and matter/hidden sector interactions (related 
to the T term) could be now investigated. In the spirit of AdS/CFT correspondence, the 
CFT generalization amounts to introducing non trivial dynamics in the bulk and brane-bulk 
energy exchange (and bulk self-interaction) in the 7D picture. The bulk cosmology with 
parameter turned on has been analyzed in section IV. 



B. Critical points and stability 



The fixed points can be derived as we have done in appendix B for the conformal non 
interacting theory and their stability can then be studied for specific theories. We won't 
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discuss this topic here. Since the new deformation parameters X, T depend on the 6D space- 
time curvature, they contain functions of the Hubble parameters and spatial curvatures and 
of the intrinsic energy scale of the background, the AdS^ radius (or Mpi). They will thus 
in general modify the equations for the fixed points and their stability in a sensible way, 
depending on the specific generalization one wants to consider. 

We will instead try to understand how the comparison with the bulk description gets 
changed when we go to the generalized scenario. This will be the subject of the next 
subsection. 



C. Comparison to 7D cosmology with energy exchange in slowly scaling regime 

As we have done for T = X = 0, wc want to illustrate some explicit examples with 
the aim of understanding the peculiar features of this cosmological model and its two dual 
descriptions. We will therefore make some assumptions simplifying the set of equations 
including Einstein, conservation and anomaly equations. First of all, we are going to neglect 
terms containing higher orders in the background curvature — namely the anomaly coming 
from Srs and the trace contribution Y coming from the second order action Sji^ . 

Equal scale factors The correspondence works as in the conformal non interacting anal- 
ysis of subsection VlllAl. On the brane side (referring to eqs (185)~(188)), the slowly 
scaling approximation leads to the equations 

X + GHx ^ HX + T (190) 

p + {3{1 + w) + 2{1 + w^))H p ^ -T (191) 

To get the bulk description expressions for H, p, x we truncate equations (35)-(37) to first 
order in the density, neglecting p/V w.r.t. order 1 terms 



X + 6i/x = - ^i/T/ (193) 

p+ (3(l + «;) + 2(l + w^))i/p = -2r° (194) 

The matching with the system of equations on the brane (189)-(191) is exact if we have 
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the following relations among the brane and bulk parameters 

T = 27;° (196) 

In the previous equations we have also explicitely evaluated the matching for equal matter 
pressures w„ = w, that gives the RS fine-tuning X^s = 0, as in the non interacting conformal 
theory. 

Static compact extra dimensions The condition of static internal space F — Q, toghether 
with the small Hubble parameter approximation, brings equations (185)-(188), relative to 
the brane description, and equations (35)-(37), relative to the bulk description, in the form 

X + ^Hx^HX + T (199) 
p + 3{1 + w)Hp = -T (200) 

(where we assumed a; = 1/3) and 

40 

X + AHx^ 2T0 - A^HT^ (202) 
p + 3(w + l)Hp = -2T^ (203) 

The parameters in the gauge and gravity descriptions are thus related by the following 
expressions 

T = 2T^ (205) 

For = (w + 5)/6 we get the zero effective cosmological constant on the RS brane, as 
before. 

Proportional Hubble parameters In the limit of proportional Hubble parameters or equiv- 
alently scale factors related by 6 = a^'', we use the set of equations (185)-(188), substituing 
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F — ^bH, K, ^ Kb, and expanding in the slowly scaling approximation 

X + = (3u; + ^6(1 - 3a;)) M + T (208) 

p + = -T (209) 

As before = 3(1 + uj) + 3^;,(1 — uj), w^^ = 3(1 + w) + 2^ft(l + w^). The bulk dynamics is 
described by (38)-(40) with F = ^bH, k ^ kb 

40 

X + d^^Ex = 2r° - -^//T/ (211) 
p-rW^^Ep = -2T^ (212) 

with = 6(e| + 2^5 + 2)7(2^5 + 3), w^^ = 3(1 + w) + 2^5(1 + w^). 

If w and -u^TT ai'c the same on both sides of the duality, then we must make the identification 
— = to match p equations. As a consequence, uj = 1/(2^ + 3) in order to have 

agreement for the mirage density (non) conservation equations and Kb = (^^+6^+3)kb/(^^+ 

3^ + 6). With these conditions, the comparison between the two sets of equations thus gives 

the following matching relations 

- ^^ + ^ (213) 

T = 2r° (214) 
2^ + 3 , X e + Q^ + ^ Tj 

2^^ + 3 2cv,iV ^ Mj>i 2^ + 3 M^o ^ ' 

X. SUMMARY AND CONCLUSIONS 

In the context of holographic cosmology, we have investigated the specific background 
of 7D RS gravity, including an energy exchange interaction between brane and bulk. Some 
novel features arise both on the bulk side of the duality and in the conformal holographic 
theory on the brane. In particular, we found distinctive results in the comparison between 
the two descriptions that need a better understanding. The originality with respect to the 
5D/4D holographic cosmology [10, 43] is due to the compactification over a 2D internal 
space, around which we wrap the 5-brane. The 6D space-time filled by the brane acquires 



62 



a dishomogeneity that distinguishes the 3D visible space from the 2D internal directions. 
Evolution can generally be different in the two spaces and pressures are individually related 
to the energy density by the usual ansatz p = wp and tt = w^p {p and tt are respectively 
the 3D and 2D pressures), with w^, w in general. 

Concerning the gravity theory in the bulk, we studied the Friedmann like equation that 
comes along with the introduction of a mirage energy density satisfying to the non homo- 
geneous radiation equation in some effective number of dimensions (which is six for equal 
scale factors in both 3D and 2D spaces and is four when the internal space is static). The 
bulk cosmological evolution is then determined by the Friedmann equation and by the (non) 
conservation equations for the mirage density and the localized matter density on the 5- 
brane. Making use of some simple ansatz for the evolution of the 2D compactification space 
(such as putting the corresponding Hubble parameter F equal to the Hubble parameter of 
the visible space H or to make it vanish) we found a wide spectrum of possible cosmologies 
that reduce to the RS vacuum in the absence of matter (i.e. we imposed the RS fine-tuning 
Xrs = 0). 

Assuming small density approximation, we have described the explicit analytical solution 
in case of radiation dominated universe. The Hubble parameter evolves as in an effective 
6D (4D) radiation dominated era for equal scale factors (static compact extra dimensions), 
independently of the form of the brane-bulk energy exchange. The effective 4D mirage 
and matter energy densities obey to the 4D free radiation equation in absence of energy 
exchange. If energy flows from the brane into the bulk, the 4D localized energy density is 
suppressed in time, in favor of the mirage density, even with zero mirage initial condition. 
For influx to the bulk, the 4D matter energy density apparently grows unbounded (if T 
is linear in p, otherwise energy density may go to zero for suitable values of the theory 
parameters), eventually diverging at a flnite time. The small density approximation must 
break down and the full analysis is needed. On the other hand, still for small densities 
but generic perfect fluid equation of state (non necessarily pure vacuum energy) and energy 
influx, we found inflationary flxed point solutions that are stable for a large class of energy 
exchange parametrizations T = Ap'^ . These thus represent stable de Sitter solutions for (the 
four visible space-time directions of) our universe. We moreover argued that, differently 
than in the 5D RS approach [43] , we may have a stable de Sitter critical point solution even 
for energy outflowing to the bulk, in the case of equal scale factors with w < —1/3 (and 
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v — 1). For dynamical compactification (i.e. F — ^H, ^ < 0) [81] we could in principle also 
get an outflow stable inflationary fixed point. We note that the 4D mirage energy density 
evolution without energy exchange is governed by the effective 4D free radiation equation 
only in the two limits of equal scale factors F = H and static compact extra dimensions 
F = 0. 

Having dropped the small density approximation, more elaborate models of cosmologies 
were developed. The number of possible inflationary critical point solutions can be larger 
than one, depending on the parametrization for the brane-bulk energy exchange. For energy 
influx we showed the 6D picture of a scenario with two fixed points, where trajectories in 
the phase space can either always be characterized by positive acceleration, either remain 
at all time with negative acceleration, or alternate acceleration and deceleration phases. 
The portraits are rigorously valid for the effective 4D energy density in the case of static 
internal space (up to a constant rescaling). If we have equal scale factors, the evolution 
equations become much more complicated functions of the 4D densities and the computation 
is beyond the scope of the paper. For u — 1 there seems to exist only the trivial critical 
point characterized by vanishing Hubble parameter, so that the energy density should grow 
without bounds as predicted by the small density approximation, until the full analysis 
is needed. For energy outflow, the 6D energy density localized on the 5-brane decreases 
and the trajectories in the phase space go toward the trivial fixed point, eventually passing 
through an accelerated era. The effective 4D picture may differ from this description in the 
equal scale factor case, since it would be possible in principle not to have decreasing density 
at all times. 

The rather detailed study of the various cosmologies emerging from the 7D RS model with 
energy exchange has been embodied in the context of the AdS/CFT correspondence. We 
have examined the role played in the holographic critical point analysis by the 6D anomalous 
CFT coupled to 6D gravity — with the addition of the higher order counterterms to the 
dual action. The dual theory on the brane is conformal (classically) and non interacting 
(with the matter theory on the brane) . This CFT would then correspond to the RS set-up 
with no energy exchange. Despite this fact, we may find infiationary critical point solutions, 
depending on the anomaly parameters and on the coefficient of the second order counterterm 
(in the curvature). All this fixed points are characterized by zero matter energy density. 
Clearly, neglecting all the higher order contributions (including the anomaly .4(6) and the 
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trace term Y) we recover the trivial fixed points of the pure RS gravity background. 

The comparison between the two dual descriptions has been achieved in the approxima- 
tion where all the higher order terms can be neglected, i.e. for small Hubble parameters. 
Since higher order terms are truncated, we cannot access to the typical non conventional 
dependence in the expression for — only linear terms are present in this approximation. 
Comparing the bulk Friedmann equation with the corresponding equation derived in the 
holographic description, we have to match the ratio M^^/V (M is the 7D Planck Mass and 
V is the tension of the RS brane) in terms of the 6D Planck mass Mpi (the 4D Planck mass 
M(4) is related to Mpi by Mf^-^ = 1^(2) M|;). The RS fine-tuning M^^ = 20M^° /V is restored 
when we recover homogeneity in the background, imposing F = H and w^^ = w. With these 
assumptions indeed, the matching is exact also with respect to the spatial curvature terms. 
As we move away from homogeneity, we have to define an effective spatial curvature for the 
compact extra dimensions in the holographic description. The matching between the scales 
of the theories refiects the RS fine-tuning for a specific value of the matter pressure in the 
internal space (determined by Wt^), depending on the proportionality factor ^ relating the 
two Hubble parameters F{t) = ^H{t), or b{t) = a^{t). We finally matched the evolution 
equations in the generalized holographic dual theory with the general bulk description. The 
interactions between hidden and visible sectors encode the dynamics of the brane-bulk en- 
ergy exchange, Ty ^ 0, on the bulk gravity side, while the breaking of conformal invariance 
(via non zero /3-functions or masses) amounts to turning on the bulk "self-interaction", 
T/ ^ 0. 

The 7D RS background has been quite accurately studied on the bulk side, though many 
profound cosmological aspects have not been explored. There could be space, however, to 
fit the cosmological evolution of the universe in this model, since one of the stable de Sitter 
critical point solutions we found could represent the actual accelerated era. Besides, trajec- 
tories can end into the stable point first passing through a decelerated phase representing the 
matter or radiation dominated universe. Another accelerated era may be present at early 
times, eventually corresponding primordial infiation. Still, there is no rigorous construc- 
tion of such a precise evolution. The holographic dual theory could also give an interesting 
cosmological description of the brane-world. It would be interesting to exploit Starobinky 
argument of graceful exit from primordial infiation via higher derivatives term in this con- 
text. In conclusion, we showed the basics of the brane-bulk duality in 7D RS background 
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with energy exchange and of its cosmological features, but many further questions can be 
adressed within the framework of the 7D RS holographic cosmology. 
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APPENDIX A: CONFORMAL ANOMALY AND TRACES IN SIX DIMEN- 
SIONS 

Conformal anomaly The conformal anomaly for 6D theories has been studied in [78]. 
It can be derived using AdS/CFT and the gravitational renormalization procedure as in 
[72, 73, 76]. 

In our notations, the general expression for the trace anomaly in a six dimensional CFT 

is 

^(6) = - (cAEiQ) + Cb/(6) + J('^5)) (Al) 

£^(6) is the Euler density in six dimensions (type A anomaly), /(g) is a fixed linear combi- 
nation of three independent Wcyl invariants of dimension six (type B anomaly) and ViJ^'^-j 
is an linear combinations of the Weyl variation of six independent local functionals (type 
D anomaly), so that at the end we have eight free coefficients in the general form of the 
anomaly, depending on the specific CFT. The type D anomaly is a trivial (it is a total 
derivative, indeed) scheme dependent term that can be cancelled by adding local covariant 
counteterms to the action [75] . 
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For our metric we obtain as a result that £^(6) depends on the Hubble parameters and on 
their time derivatives up to order one (that is, up to the second time derivative of the scale 
factors). J(6) instead depends on H and F time derivatives up to order three, and so does 
the divergence term. To be more specific, /(g) is made up by three contributions /i,/2,/3, 
with fixed coefficients; the first two are two different contractions of three Weyl tensors (and 
contain only derivatives of the Hubble parameters up to order one), while in /a there are 
second order derivatives of the Weyl tensor (i.e. third order time derivatives of the Hubble 
parameters) . 

For a 6D FRW background (i.e. requiring homogeneity in all six dimensions) the sum of 

type A plus type B anomalies depends on the Hubble parameters only up to the first time 
derivative, while the type D anomaly contains time derivatives up to order three. 

In terms of the Riemann tensor, Ricci tensor and scalar curvature, the A, B and D 
contributions to the anomaly read [72] 

^(6) = 

1/10. 1 1 



^^^^ 1152 V 3 6^^ ^ 10^^ 

^ -R^"'V^R + -^RV^R (A2) 



2880 5760 



where 



Eo = K: 



- I2K2 + 3K3 + I6K4 - 24X5 - 24:Kq + 4i^7 + 8Ks 

q 27 3 5 3 

^ 200 40 10 ^ 4 ^ 2 ^ ^ ^ 

h^Ki- SK2 - 2K3 + 10X4 - 10X5 - ^^9 + 5X10 - 5X11 



and 



ixpcrv 1 



(Xi, . . . , Kii) — (i?^, RR^^R'^", RRn^pfrR'^"'"^ , R^^ Ru'' Rp'^ ^ R'^^R'^'^R^ 

-"'^i'-"' -'I' pari ^^p,v per rm ■'■^pupa^^ tA ) 

RUR, Rp^nRf^", Rp.p.UR^^'^-) 

In the analysis of the solutions to the Friedmann equations we plug in the specific expression 
for the Riemann tensor obtained considering our ansatz (108) for the metric. But before 
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doing this, we use the anomaly equation and some standard assumptions on the pressures 
that parametrize the stress-energy tensors to manipulate our system of differential equations. 

To give an explicit result for the conformal anomaly in the specific case of 6D CFT on 
curved space-time, with the ansatz (108) for a 4D FRW plus a 2D compact internal space 
background, we write the type A contribution, in terms of the 3D and 2D spaces Hubble 
parameters H, F and spatial curvatures k, k: 

+2 (f + F^) (h^ + ^) } (A3) 

The type B and D contributions have a more complicated form and we write them when it 
is necessary, in the specific limits we consider throughout the paper. 

For the (0,2) SCFT dual to the M5 background, the anomaly coefficients are given by 
CA^CB^myn^ [72]. 

Counteterm traces The dual RS theory action contains the three counterterms Si, S2, 
S3 written in (96). Varying these contributions w.r.t. the six dimensional induced metric 
7^,^ on the brane yields [98] [76] 



1 

32 



2 3 
5 5 



-^7,. {RpaR'" - - l^R^ - t;, log ^ 



(A4) 



Where T^^ is a traceless tensor of cubic order in the curvature. The trace of the conformal 
variation of 5*1 (corresponding to the linear part of (A4) in the curvature) gives a term 



proportional to the Einstein tensor. The variation of the S3 (x J ^(0)0(6) action (related 
to T^^, where we introduced a(3) as in the standard notation of [73]) is traceless because it 
is proportional [76] to the traceless tensor /i(6);ui/ that enters into the parametrization of the 
metric (91) due to Fefferman and Graham. However, the variation under conformal trans- 
formations of the cutoff dependent counterterm is non trivial and is the only contribution 
to the conformal anomaly, so that A(^q) oc a(g) (see for instance [76, 91] for a more detailed 
derivation). Finally, the trace of 5*2 (equal to the trace of the quadratic contributions in 
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(A4)) is 



Y = —MH ( R'^'Rau - — 

32 V ^ 10 J 



(A5) 



which can be expressed in terms of the Hubble parameters H, F and of the spatial curvatures 

k, K of the (4D FRW + 2D compact space) background (108) as 
2M^i 



Y 



160 



3^ + 64 (3^ + F'' + 8FH + 3H' + 6F + 4h') + 



a2 V 62 

(F - 6H) {F + 3H) + F + 9H 



+ 



-3F^ + ASF^H + 6FH (21H'^ + 7F + 13H 



niH^-AF + 2AH^ + 



+3 



6H' 



- (^F-H^ + 2H^ 



+ 2H^ 17F + 3H 



')]} 



(A6) 



We define cy = M^i/32Mpi wich is given as a function of the number N of M5-branes 
in the gravity background by cy = \/2N^Jt^. 



APPENDIX B: FIXED POINTS IN THE HOLOGRAPHIC DESCRIPTION 

In this appendix, we are going to look for the existence of inflationary points for our 
specific holographic model and to find what kind of dependence they have on the parameters 
of the theory. We will also study the stability matrix determining — in some special limits 
— whether the critical points are stable or saddles. 

In the calculations, we suppose that the effective cosmological constant on the brane A is 
null. 



1. Flat compact extra dimensions 

We start by considering the case of (locally) flat internal space, which could be for example 
a two-torus. The three spatial dimensions of the 4D FRW are already supposed to be flat, 
so that the system of equations of motion simplify, having dropped the terms proportional 
to both spatial curvatures. 

The general flat extra dimension fixed points (F* ^ 0) are not easy to characterize. 
We choose to analyze the case in which the extra dimensions Hubble parameter is zero at 
the fixed point, meaning that the fixed point represents a universe with static flat extra 
dimensions. 
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a. F^^O,uj^l 



Fixed point solutions Looking for the solution to the Priedmann plus conservation set 
of equations with constant Hubble parameters {H = H^, F = F^) and zero curvatures 
{k — K — 0), we have to consider the simplified system of equations (where we have already 
solved the equation for a) 

{3H'^ + 6H, F, + F^) - [SuH, + (1 - 3a;) (^(e). + F.) = A (Bl) 

{3Hl + AH, F, + 3F2) - u + (^(g), + n) = A (B2) 

a, = [SuH, + (1 - So;) F,] A^e). = {SH^ + 6H,F, + F^) - A (B3) 

= -u + 2F,) (i(6)* + n) = -M^, {SH^ + AH^F, + 3F^) + A (B4) 

p, = 0, x* = (B5) 

where the relation between Q and a; was defined to be 1/a; = 2^2 + 3 and the trace contri- 
butions ^(6)*, y* take the form 

+291F^Hl + 18AF^H!^ + bb7F^H^ + 138FH^ - MH^] (B6) 
6 

n = ^cyM^pi {2Hl + 2H,F, + F^) {3Hl + 18H,F, - F^) (B7) 

and ^(6)*, have been defined as 

i(6)* = ^(6)*/ [3(1 + uj)H, + 3(1 - uj)F,] 
n = n/ [3(1 + u)H, + 3(1 - u)F,] 

for (1 + u;)//^ + (1 - uj)F, ^ 0. 

For F, y^O and u; 1/5, we can reformulate eqs (Bl) and (B2) in order to get 

(3 - 21u;)Hl + (4 - 30)F,H, - (3 - l^F^ = (1 - 5u) (B8) 

(^(6)* + n) = 2M|, (B9) 

Imposing A = (no effective constant on the brane) greatly simplifies the solution since 
H^, oc F^,. Under that assumption, defining Ce and — as functions of ui and e = ±1 (Pg 
is a function of the anomaly parameters ca and cb, of cy and the Planck mass as well) — 
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such that — — CeF and ^«4(6)* + j = T^eF^, the solution takes the form 



H'. = Ml, (C, + if 



2a 



(1 - buj)V^ 



, Ft = Mil 



2a 



(BIO) 



This solution exists for the values of lo such that Ce/V^ > (for u; < 1/5) or C^/V^ < (for 
CO > 1/5). 

The CFT energy density and pressures are then given by 



2u; 



2a 



(1 - 5uj)V, 



(for uj = 1/3 we have a^^ = 0). 



= -M^;a;(5 + 3a) T), 



2Ce 



(1 - 5u;)V, 



(Bll) 
(B12) 



Fixed point solutions To analyze the case cu — 1/5, it's better to reformulate equations 
(Bl) and (B2) in the following way: 



2M; 



PI (H^ - F^) F^-(l- hu)F^ (i(6H + n) = (B13) 

{ml + F, + Fl) - [3ujH, + (1 - Soo) F,] (^(6). + n) = A (B14) 

Prom the first equation we get — F^ and substituing in the second we find the equation 
for //^ 

5 



288 



[ca + '2cb)H° - UOcyMliH^ + lOM^^//^ = A 



(B15) 



For A = it has a non trivial solution only if (24cy)^ > {ca + ^cb) 

24 



Ht = Ft 



-Mil 



Ca + 2cb 

The energy density and pressures are equal to 

240 



24cy ± a/5764 + (ca + 2cb) 



(B16) 



-a. 



p* 



Ca + 2cb 



Mil 



24cy ± ^5764 + {ca + 2cb) 



(B17) 
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c. F^ = 



Fixed point solutions Supposing instead F^, — 0, the fixed point solution is 



20 uj 



3cb lj + 1 

2lo 



Ml, 



48cy ± ^/6 ( 3844-cb 



UJ 



(7+ = -a. 



3a; — 1 



20 u 
cbu; + 1 



48cy ±4/6 3844- cb 



cu+l 



(B18) 

(B19) 



for u ^ —1. This gives real Hubble parameter for 384cy — cbuj/{}jJ + 1) >0. Ifa; = —1, we 
find 



640c 



Y 



640c, 

cb ' " '2 3cb 



Mil, 







(B20) 



If the CFT is characterized by a positive Cb, there is no non trivial critical point. For 
vanishing cb the only fixed point with A 7^ is the trivial one. 

When F^ — —H{1 + a;)/(l — a;), there is only one possible solution, for which the param- 
eters must have the values: cu — —1 (i.e. a^^ — —2ap), cb — 0, F^, — and the fixed point is 
thus the one in (B20). 

Stability analysis For both fixed points characterized by the zero value of the extra 
dimension Hubble parameter, i.e. both for co ^ —1 ot 00 = —1, we must find the eigenvalues 
of a 4 X 4 matrix. In fact we have a third order linearized difi'erential equation for the 
perturbation SH{t) and a first order ODE for the energy density dp, while SF{t) is found to 
be proportional to SH{t) solving an algebraic equation 



SH^^'^ = -a2 5H -ai5H- qq SH + cq Sp 
Sp = -3{l+w)HJp 
SF^aSH 



(B21) 
(B22) 

(B23) 



The coefficients in the differential equations are functions of the anomaly parameters ca,cb, 
of the trace parameter cy, of a; and of Mpi. 

The eigenvalues Ai, A2, A3 are then given by the roots of the third degree polynomial 



A^ + a2A^ + aiA + ao = 



(B24) 
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while A4 = —3(1 + w)H^ < 0. The coefficients oq, oi, 02 are given by 

12 8000M^,(3 + 2q; + 3(1 + aW) , ,,4,^^ , .2 r.2N 

1 960000M^,(H-a) , ......^ ^^^v , ,2 ..2n 

= ^^i/. (B25) 

1 — a 

where 

800M|,(l + c^)-(-3c^if,^ + 480cyM|,ff,^)o. 
800M^,(5 - 3c.') - (-3cbH^ + A&QcyMIiHI) (1 - 3cu) ^ 

The sign of the eigenvalues Ai, A2, A3 (A3 = A2 \S.27a'l + Aa\ — lSaoaia,2 — a\a2 + AaQa\ > 0, 

otherwise we get three real roots) determines the nature of the fixed point. Since A4 < 0, 

we find that we can only have a completely stable fixed point or a saddle. In the case of one 

real and two complex conjugated roots the critical point can be attractive only if 02 > and 

^^^<A + 5<^ (B27) 



3 3 

where 



A - sgn(i?) (\R\ + ^R^ - g3 j \ B = ^ (B28) 

R=^ {2al - 90102 + 27oo) , Q = ^ (o^ - 3oi) (B29) 
54 9 

When the three roots are real, they are negative (corresponding to an attractive fixed point) 
iff Oqi Q-ii Q'2 > 0. For the other values of Oq) Q-i, ^2 the critical point is a saddle. 

The coefficients of the linearized differential equation don't depend on the anomaly pa- 
rameter ca corresponding to the type A anomaly, so that only type B anomaly infiuences 
the characteristics of this fixed point. 



2. Static compact extra dimensions 

We analyze the set of differential equations when the extra dimensions are (locally) com- 
pactified on a sphere, i.e. k > 0, supposing that the corresponding acceleration factor b{t) 
remains constant, so that F(t) = 0. 

Beside the — fixed points, we only have two acceptable time independent solutions 
to the Priedmann equations. The — 0, k ^ fixed points are always saddle points as we 
can conclude from the linear order analysis, since the eigenvalues of the stability matrix (or 
their real parts) are one opposite to the other. 
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a. u ^ 



Fixed point solutions The energy density is zero, if* and cr* are then determined by 

Cb 



4800 



bo bo 



+ 



4c 



PI 



M 



PI 



m + 12 



K ^ 







+ 



(B30) 

(B31) 
(B32) 



Restricting the possible values of a;, we can obtain at least one positive root of eq (B30), 
without having limitations on the anomaly coefficients ca,Cb- 

We can illustrate an example, choosing the simple case cb = (i.e. there is no contribu- 
tion from the conformal invariants in the anomaly) and also u — —1, which simplifies the 
equation (B30). The fixed point is thus determined by 

-1 



9± 4 90 



5ca 1 K, 
192cy Mj,i 6| 



= 



2a. 



Mj>i I 3 



9± W90 



5ca 



K, 



192cy bl 




(B33) 
(B34) 



which is real for 192cyM|,; > — cak/456q. We can moreover have two distinct positive 
fixed points if 5cA/t/96g < -192cf.M|,;. 

Stability analysis We now analyze the H^, ^ fixed points behavior. 

Regarding the fixed point determined by (B30), we get a negative eigenvalue A4 = —3(1 + 
w)H^ (given that w > — > 0) and the other three are the roots of the third degree 
polynomial 



+ a2X'^ + aiX + ao = 



(B35) 
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where = ai/a^, i — 0,1,2 and 

We get a 4 X 4 stability matrix — despite the fact that we should have only 2 variables 
{H, p) — because the differential equations are of third order: the p eigenvalue is A4, but H 
is a superposition of the four modes corresponding to the four eigevalues of the matrix. 

As in the previous analysis for the flat extra dimensions, the solutions Ai^2,3 of the equation 
(B35) are such that Ai e M, A3 = A2 or Ai, A2, A3 e M. Besides, when we have the complex 
conjugated pair, there are only three possibilities: 

1. Ai, 3?(A2) = 3?(A3), A4 < =^ the solution is stable (even if one of the eigenvalues are 
null, because that mode won't then contribute to the expression for H) 

2. Ai, A4 < 0, 3?(A2) = 3?(A3) > ^ we get a saddle point 

3. Ai > 0, 3?(A2) = 3?(A3), A4 < in this case too, the fixed point is a saddle 

The equahties Ai = and 3?(A2) = 3?(A3) in the case (i) are possible, but not simultaneousely. 
When the roots are all real we can get a stable point iff ao, oi, 02 > 0, which implies 
Ai, A2, A3 < 0, and a saddle otherwise, with one negative two positive, or two negative one 
positive roots. 

3. Equal scale factors 

Another limit that simplifies some of the calculations is the equal scale factor assumption. 
In this case the Hubble parameters of the internal space and the 3D space are equal, F — 
H, and we remain with a set of equations for the variables H, p,a, as in the static extra 
dimension limit. 
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a. u ^ 



Fixed point solutions We observe that when u — there is no acceptable solutions to 
the time independent Einstein equations. So, we want to find the fixed points in the H — F, 
cu ^ limit. The time independent Priedmann plus conservation equations lead to 



6M 



(l0//,2 + ^) (B37) 

a, = Mj,i {lOHl + ^) (B38) 
= (B39) 

provided that 5 + 3w + 2^^ 7^ [99] . As a consequence of solving the system of equations 
wc also obtain a constraint on the CFT pressure of the hidden sector cTj^, as wc must impose 
u = 1/5, i.e. (7^ = Cp. The equation (B38) yields the value of the Hubble parameter at the 
fixed point as a function of ca,cb,cy, Mpi, k. 

In the case of fiat extra dimensions k = we immidiately solve the system of equations 
finding (discarding the trivial — solution) 

24cy ± ^5764+ (cA + 2cs) 

24cy d= ^5764+ (ca + 2cb) 

The solution is acceptable if —{ca — 2cb) < (24cy)^ and the two roots are both positive 
when {ca — cb) < 0. 

Stability analysis The last situation that we are considering is the case of equal scale 
factors, as in the fixed point analysis. We only found one fixed point with F = H, that 
entails a relation between the two CFT pressures a-j^ = o'p {u = 1/5). We could thus calculate 
the stability matrix eigenvalues corresponding to this particular limit. 

When the extra dimensions spatial curvature is zero k — 0,m addition to the vanishing 3D 
curvature {k — 0), the stability matrix can be studied straightforward. All the eigenvalues 
are coincident since 5H oc 5a oc 5p. They are given by 

X = -{5 + 3w + 2w^)H^<0 (B41) 

For if* > the fixed point is hence stable. 



24 

CA + 2CB ^' _ 

240 

c^* = -CTp* = -cTtt* = —: — 

Ca + ^cb 



(B40) 
P* = 
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